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Abstract 

In this review we consider two different models of a hydrogenic atom 
in a quantized electromagnetic field that treat the electron relativistically. 
The first one is a no-pair model in the free picture, the second one is 
given by the semi-relativistic Pauli-Fierz Hamiltonian. For both models 
we discuss the semi-boundedness of the Hamiltonian, the strict positiv- 
ity of the ionization energy, and the exponential localization in position 
space of spectral subspaces corresponding to energies below the ioniza- 
tion threshold. Moreover, we prove the existence of degenerate ground 
state eigenvalues at the bottom of the spectrum of the Hamiltonian in 
both models. All these results hold true, for arbitrary values of the fine- 
structure constant, e 2 , and the ultra-violet cut-off, and for a general class 
of electrostatic potentials including the Coulomb potential with nuclear 
charges less than (sometimes including) the critical charges without ra- 
diation field, namely e~ 2 2/TY for the semi-relativistic Pauli-Fierz operator 
and e _2 2/(2/7r + 7r/2) for the no-pair operator. Apart from a detailed dis- 
cussion of diamagnetic inequalities in QED (which are applied to study 
the semi-boundedness) all results stem from earlier articles written by the 
authors. While a few proofs are merely sketched, we streamline earlier 
proofs or present alternative arguments at many places. 
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1 Introduction 

In the late 90 's and the past decade the existence of ground states of atoms 
and molecules interacting with the quantized photon field has been intensively 
studied by mathematicians in the framework of non-relativistic quantum electro- 
dynamics (QED). The corresponding Hamiltonian is the non-relativistic Pauli- 
Fierz operator which, in the case of a hydrogen-like atom, is given as 

H™ := (a ■ (-iV x + A)) 2 - SL. + H f . (1.1) 
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Here <r = (oi, 02, 03) is a vector containing the Pauli spin matrices, A is the 
quantized vector potential in the Coulomb gauge, and Hf is the energy of the 
photon field. The symbol A includes a prefactor entering into the analysis as a 
parameter, namely the square-root of the fine structure constant which equals 
the elementary charge, e > 0, in the units chosen in this paper. The Coulomb 
coupling constant, 7 > 0, is the product of e and the nuclear charge. We shall, 
however, always consider it as an independent parameter since the interrela- 
tionship between e and 7 does not play any role in our work. A additionally 
depends on some ultra-violet cut-off parameter, A > 0. By now it is well-known 
that has a self-adjoint realization in the Hilbert space L 2 (R"^, C 2 ) <g> J^bP^] 
whose spectrum is bounded below. Here \%\ denotes the bosonic Fock space 
modeled over the Hilbert space for a single photon, Jif — L 2 (R 3 x Z2). Proving 
the existence of ground states for H™ T means to show that the infimum of the 
spectrum of iJ" r is an eigenvalue corresponding to some normalizable ground 
state eigenvector in £ 2 (R3, C 2 )®^b[«#1- Because of the spin degrees of freedom 
this ground state eigenvalue will be degenerate. Mathematically, the study of 
the eigenvalue at the bottom of the spectrum of H^ T is very subtle because the 
spectrum of H™ is continuous up to its minimum and the eigenvalue is, thus, 
not an isolated one. In particular, many standard methods of spectral theory 
do not apply and several new mathematical techniques had to be invented in 
order to overcome this problem. 

The first proofs of the existence of ground states for H" r and its molecular 
analogs have been given in [4j [6] , for small values of the involved parameters e 
and A. A few years later the existence of ground states for a molecular non- 
relativistic Pauli-Fierz Hamiltonian has been established, for arbitrary values of 
e and A, in }17j by means of a certain binding condition which has been verified 
later on in 33 . Moreover, infra-red finite algorithms and spectral theoretic 
renormalization group methods have been applied to various models of non- 
relativistic QED to study their ground state energies and projections [HEHHO 
O [7J [TS] . These sophisticated methods yield very precise results as they rely on 
constructive algorithms rather than on compactness arguments as in [?J [SJ [T7J . 
They work, however, only in a regime where e and A are sufficiently small. 

A question which arises naturally in this context is whether these results 
still hold true when the electrons are described by a relativistic operator. In 
this review, which summarizes results from [25l [28l EH EH |37j [48] , we give a 
positive answer to this question. We study two different models that seem to 
be natural candidates for a mathematical analysis: The first one is given by the 
following no-pair operator, 

:=P+(P A -^ + P f )P+, (1.2) 

or more generally, 

Hy P ■= Pa {D a + V + H f ) P£ , (1.3) 

for some electrostatic potential V. Here Da is the free Dirac operator minimally 
coupled to A and denotes the spectral projection onto the positive spectral 
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subspace of Da, 

■f a : = l(o,oo) (Da)- 
The no-pair operator is considered as an operator acting in the projected Hilbert 
space, J4?j£ := Pt Stf, It is thus acting on a space where the electron and 
photon degrees of freedom are always linked together. The analog of H^ p for 
molecules has been introduced in [32] as a mathematical model to study the 
stability of matter in relativistic QED. Under certain restrictions on e, A, and 
the nuclear charges it has been shown in [3 2) that the quadratic form of a 
molecular no-pair operator is bounded from below by some constant which is 
proportional to the number of involved electrons and nuclei and uniform in the 
nuclear positions. Moreover, the (positive) binding energy has been estimated 
from above in [31 . In fact, there are numerous mathematical contributions 
on no-pair models where magnetic fields are not taken into account or treated 
classically; see, e.g., [39 for a list of references. For instance, it is shown in [13] 
that a no-pair operator with Coulomb potential but without quantized photon 
field - which is then often called the Brown-Ravenhall operator - has a critical 
coupling constant, 7" p := 2/(2/7r + 7r/2), such that the corresponding quadratic 
form becomes unbounded below when 7 exceeds this value. Moreover, various 
molecular no-pair models (without quantized fields, however) are widely used in 
quantum chemistry and in the theoretical and numerical study of highly ionized 
heavy atoms; see, e.g., [THIH^, 46 . In this context several different choices of the 
projections determining the model find their applications. For instance, one can 
include the Coulomb or a Hartree-Fock potential in the projection. These two 
choices are covered by the results of [39] where two of the present authors provide 
a spectral analysis of a class of molecular no-pair Hamiltonians with classical 
magnetic fields. The choice of the projection P£ which does not contain any 
potential terms is referred to as the free picture. We remark that it is essential 
to include the vector potential in the projection determining the no-pair model. 
For, if Pj£ is replaced by Pq~, then the analog of (|1.2[) describing N interacting 
electrons becomes unstable as soon as N ^ 2 [TB l l32 ] l3"3]. Moreover, the operator 
in (ll.2[) is formally gauge invariant and this would not hold true anymore with 
Pq in place of P^. Gauge invariance plays, however, an important role in the 
proof of the existence of ground states as it permits to derive bounds on the 
number of soft photons. In fact, employing a mild infra-red regularization it is 
possible to prove the existence of ground states for the operator in (|1.2j) with 
P£ replaced by P + [Ml [35]. It seems, however, unlikely that the infra-red 
regularization can be dropped in this case [24] . 

The second operator treated in this review, the semi-relativistic Pauli-Fierz 
operator, is given as 

v/(<r-HV + A)) 2 + l-^+ff f , (1.4) 

where <x is a vector containing the Pauli spin matrices. Since \f— A and l/|x| 
both scale as one over the length there will again be some critical upper bound 
on all values of 7 > for which (|1.4I) defines a semi-bounded quadratic form. 
As we shall see this upper bound is at least as big as (in fact equal to [28]) 
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the critical constant in Kato's inequality, 7 PF := 2/ir. Again we shall study 
the semi-relativistic Pauli-Fierz operator also for a more general class of elec- 
trostatic potentials V. The latter (straightforward) generalization is relevant in 
a forthcoming work of the first two authors devoted to the enhanced binding 
effect [37]. 

Also the semi-relativistic Pauli-Fierz operator has been investigated earlier in 
a few mathematical articles. For instance it appears in the mathematical study 
of Rayleigh scattering [TJ] where the finite propagation speed of relativistic 
particles turns out to be an advantageous feature in comparison to models of 
non-relativistic QED. (The electron spin has, however, been neglected in [TI].) 
For 7 = 0, the fiber decomposition of (|1.4[) with respect to different values of 
the total momentum has been studied in [40]. Moreover, there is a remark in 
[40] relevant for us saying that every (speculative) eigenvalue of the operator in 
(|1.4p is at least doubly degenerate since the Hamiltonian commutes with some 
anti-linear involution. The existence of the renormalized electron mass in the 
semi-relativistic Pauli-Fierz model, i.e. twice continuous differentiability of the 
mass shell in balls about zero, is proved in [35], for small values of e. The 
last author has shown [48] that, when the speed of light, c, is re- introduced 
as a parameter and 7 £ [0,7 PF ), then the operator in (|1.4p converges in norm 
resolvent sense to the non-relativistic Pauli-Fierz operator in (ll.lj) . as c tends 
to infinity. Finally, there is a contribution [21] on the existence of binding in 
the semi-relativistic Pauli-Fierz model; see Remark 14.21 

We should also mention that the existence of ground states in a relativistic 
model describing both the photons and the electrons and positrons by quantized 
fields has been studied mathematically in 8 . To this end infra-red and ultra- 
violet cut-offs for the momenta of all involved particles are imposed in the 
interaction term of the Hamiltonian considered in [5] . 

In the remaining part of this introduction we explain the organization of this 
review article and summarize briefly our main results. In Section[2jwe recall the 
definitions of some operators appearing in QED and introduce the no-pair and 
semi-relativistic Pauli-Fierz Hamiltonians more precisely. Although the general 
strategy of our whole project relies on the methods developed in [1] [5] H7J the 
spectral analysis of the operators treated in this article poses a variety of new 
and non-trivial mathematical obstacles which is mainly caused by their non- 
locality. In fact, both operators do not act as partial differential operators on the 
electronic degrees of freedom anymore as it is the case in non-relativistic QED. 
In this respect the no-pair operator is harder to analyze than the semi-relativistic 
Pauli-Fierz operator since also the electrostatic potential and the radiation field 
energy become non-local due to the presence of the spectral projections Pj£ . In 
the last subsection of Section [3J we explain a few mathematical tools used to 
overcome some of the problems posed by the non-locality thus preparing the 
reader for the proofs in the succeeding sections. 

In Section [3] we provide some basic relative bounds and study the semi- 
boundedness of the Hamiltonian in both models under consideration. We start 
with a discussion of various diamagnetic inequalities for quantized vector poten- 
tials. They are employed to prove that the quadratic form of the semi-relativistic 
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Pauli-Fierz operator is semi-bounded below on some natural dense domain, for 
a suitable class of potentials including the Coulomb potential with coupling- 
constants 7 £ [0,7 PF ]. For the no-pair operator we obtain similar results with 
Coulomb coupling constants 7 £ [0, 7" p ). As a consequence, both operators can 
be realized as self-adjoint operators in a physically distinguished way by means of 
a Friedrichs extension. We point out that the results on the semi-boundedness, 
as well as all further results described below hold true, for arbitrary values of e 
and A. 

Section @] is devoted to the study of binding. For both models treated here 
we show that the infimum of the spectrum of the Hamiltonian with appropriate 
non-vanishing potential is strictly less than its ionization threshold which, by 
definition, is equal to the infimum of the spectrum of the Hamiltonian without 
electrostatic potential. To this end we employ trial functions which are ten- 
sor products of electronic and photonic wave functions and work with unitarily 
equivalent Hamiltonians in order to separate the electronic and photonic degrees 
of freedom. The unitary transformation used here represents the free Hamil- 
tonian (V = 0) as a direct integral of fiber operators with respect to different 
values of the total momentum. 

Typically, proofs of the existence of ground states in QED require some 
information on the localization of spectral subspaces corresponding to energies 
below the ionization threshold (or at least of certain approximate ground state 
eigenfunctions). Here localization is understood with respect to the electron 
coordinates. We establish this prerequisite in Section [5] by adapting some ideas 
from [4] [16]. In this section we present streamlined versions of some of our 
earlier arguments from 137| . Moreover, we implement later improvements [25j 
on parts of the results of [37] by providing optimized exponential decay rates in 
the case of the semi-relativistic Pauli-Fierz operator that reduce to the typical 
relativistic decay rates known for the electronic Dirac and square-root operators, 
when the radiation field is turned off. The class of potentials allowed for in 
Section [5] covers Coulomb potentials with 7 £ [0,7 PF ] in the case of the semi- 
relativistic Pauli-Fierz model and with 7 £ [0,7" p ) in the no-pair model. It is, 
however, possible to prove the exponential localization for the no-pair operator 
with Coulomb potential also in the critical case 7 = 7" p by another modification 
of the arguments [25] . 

The main results of our joint project are the proofs of the existence of ground 
states for the semi-relativistic Pauli-Fierz and no-pair operators. As already 
stressed above our proofs work, for arbitrary values of e and A and for a class of 
potentials including the Coulomb potential with 7 £ (0,7 PF ) and 7 £ (0,7^ p ), 
respectively. Starting from these results it is actually possible to prove the 
existence of ground states also in the critical cases 7 = 7 PF and 7 = 7" p , 
respectively, by means of an additional approximation argument. We refrain 
from explaining any details of the latter in the present article and refer the 
interested reader to [25] instead. 

The proofs of the existence of ground states given here are divided into two 
steps: 

First, one introduces a photon mass, m > 0, and shows that the resulting 
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Hamiltonians possess normalized ground state eigenf unctions, <fi m [U[5J[T7]. In 
this first step, which is presented in Section [6j we employ a discretization of 
the photon momenta as in [B]. Roughly speaking, by discretizing the photon 
momenta one may replace the Fock space J^b \<%\ by a Fock space modeled over 
some l 2 space. As a consequence the spectrum of the radiation field energy be- 
comes discrete and one can in fact argue that the total Hamiltonian has discrete 
eigenvalues at the bottom of its spectrum when all small photon momenta are 
discarded. At this point we add a new observation based on the localization 
estimates to the arguments of [5J which allows to carry through the proof, for 
all values of e and A. (In [H|B] these parameters were assumed to be sufficiently 
small.) Another technical tool turns out to be very helpful in order to compare 
discretized and non-discretized Hamiltonians (or those with and without photon 
mass), namely, certain higher order estimates allowing to control higher powers 
of the radiation field energy by corresponding powers of the resolvent of the total 
Hamiltonian. For the semi-relativistic Pauli-Fierz operator such estimates have 
been established in [TJ. In [35J one of the present authors re-proves the higher 
order estimates from |14) by means of a different and more model-independent 
method which also permits to derive higher order estimates for a (molecular) 
no-pair operator for the first time. We discuss these higher order estimates in 
Subsection l6.4l but refrain from repeating their proofs. We remark that many of 
the arguments presented in Section [5] are alternatives to those used in ■ 

The second step in the proof of the existence of ground states comprises 
of a compactness argument showing that every sequence {4> mj } with mj \ 
contains a strongly convergent subsequence. In fact, one readily verifies that 
the limit of such a subsequence is a ground state eigenfunction of the original 
Hamiltonian with massless photons. This step is performed in Section [8] in 
parts by means of arguments alternative to those in [281 |2"5] . The compactness 
argument requires, however, a number of non-trivial ingredients. First, we need 
two infra-red bounds, namely a bound on the number of photons with low 
energy in the eigenfunctions (f> m (soft photon bound) [SI E] and a certain bound 
on the weak derivatives of <j) m with respect to the photon momenta (photon 
derivative bound) [T7]- To derive the infra-red bounds one can either adapt a 
procedure proposed in [T7] (this is carried through in earlier preprint versions 
of 28. 29 available on the arXiv) or establish a formula for a(k) <p m by means 
of a virial type argument and infer the bounds from that representation. We 
outline the proof of the latter formula and of the soft photon bound for the 
semi-relativistic Pauli-Fierz operator in Section[7l The photon derivative bound 
and the infra-red bounds for the no-pair operator are derived by very similar 
procedures and we refer the interested reader to our original articles [28] [29] for 
the rather dull technical details. The arguments presented in Section [7] are also 
intended to emphasize the role of the gauge invariance of the models treated 
here. In fact, one first applies a unitary operator-valued gauge transformation 
(Pauli-Fierz transformation) and the infra-red bounds are then derived in the 
new gauge. Without the gauge transformation one would encounter infra-red 
divergent integrals. 

As soon as the infra-red estimates are established, the soft photon bound 
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and the exponential localization estimates show that the eigenvectors (f> m are 
localized uniformly in m with respect to the electron and photon coordinates and 
that their components in all but finitely many Fock space sectors are negligible. 
Moreover, the photon derivative bound implies that their weak derivatives with 
respect to the photon momenta are uniformly bounded in a suitable L p -space 
and since their energies are uniformly bounded we also know that the vectors 
have uniformly bounded half-derivatives with respect to the electron coordinates 
in L 2 . It is an idea of |17j to exploit such information by applying compact 
embedding theorems for Sobolev-type spaces to single out subsequences that 
converge strongly in L 2 . In the semi-relativistic setting considered in Section [5] 
some classical embedding theorems by Nikol'skii turn out to be useful substitutes 
for the Rellich-Kondrachov theorem employed in [17] . At the end of Section [8] 
we also show that the ground state energies of both Hamiltonians are degenerate 
eigenvalues. 

At last, in Section [§1 we present the proofs of some technical results we have 
referred to in earlier sections so that most parts of this review become essentially 
self-contained. 

2 Definition of the models 

In order to introduce the models treated in this article more precisely we first 
fix our notation and recall some standard facts. 

2.1 Operators in Fock-space 

The state space of the quantized photon field is the bosonic Fock space, 

oo 
n=0 

It is modeled over the one photon Hilbert space 

,Jt := L 2 (R 3 x Z 2 ,dk), [ dk := V / d 3 k. 

J AGZ2 

k = (k, A) denotes a tuple consisting of a photon wave vector, k £ R 3 , and 
a polarization label, A e T%. Moreover, := C and ^^[J(f] is the 

subspace of all complex-valued, square integrable functions on (R 3 x Z2)™ that 
remain invariant under permutations of the n € N wave vector/polarization 
tuples. The subspace 

% ■= C ffi C* ((R 3 x Z 2 ) n ) n J^ n) [JT] (Algebraic direct sum.) 

is dense in J?b[^H- The energy of the photon field, Hf := dT(ui), is given as the 
second quantization of the dispersion relation u>(k) := |k|, k = (k, A) £ R 3 x Z 2 . 
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We recall that the second quantization of some real-valued Borel measurable 
function w is given by (dT(zu) tpY°^ = and 

n 

(dT(w)^ n \ki, • • • ,kn) = ^^(fc J -)V' (n) (fci, ■ • • ,k n ) , neU, 

3=1 

for all V = (^ (n) )n=o e ^bW such that ([dT(n7) V] (n) )~=o G ^b[^1- We 
further recall that the creation and the annihilation operators of a photon state 
/ € J?T are given, for n € N, by 

(at(/) ^^(fcr, . . . , k n ) = n- 1 ^ f(kj) ^("- 1 >(. . . , kj-ukj+t,. . .) , 

j'=i 

(a(/)V) ( ™ _1) (fci,...,fc„-i) = n 1 / 2 Jj(k)^ n Hk,k 1 ,...,k n - 1 )dk, 

and (a f (/)?A) (0) = 0, o(/) Q = 0, where fi := (1,0,0,... ) G J^ b [^] is the 
vacuum vector. We define at(/) and a(/) on their maximal domains. For 
/, g e the following canonical commutation relations hold true on 

[a(/) , a(g)} = [at (/) , a t (a)] = , [a(f) , at (<,)] = ( / | g > 1 . 

For a three-vector of functions f = / (2) , / (3) ) € JT 3 , the symbol a"(f) 

denotes the triple of operators a J (f) := (a^/W), a"(/( 2 >), a"(/( 3 ))), where a" is 
always either a or at. 

2.2 Interaction term 

Next, we describe the interaction between four-spinors and the photon field. 
The full Hilbert space underlying our models is 

:= L 2 (R 3 ,C 4 )® J?b[~#1. 

It contains the dense subspace 

9 := C£°(R 3 , C 4 ) <g> % ■ (Algebraic tensor product.) 

We introduce the self-adjoint Dirac matrices 011,012,013, and /3 that act on the 
four spinor components of an element from Jff, that is, on the second tensor 
factor in Jf? = L 2 (R 3 ) <g> C 4 ® & h [JXr]. They are given by 

«,:=(° *'), i€{l,2,3}, ^=ao:=(; ^ , 

where cti, 02,03 denote the standard Pauli matrices, and fulfill the Clifford al- 
gebra relations 

a l a J +a j a l = 25 ij t, i,j e {0, 1, 2, 3} . (2.1) 
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The interaction between the electron/positron and photon degrees of freedom 
in the Coulomb gauge is given as a ■ A := oil A^ 1 ' + «2 A^ + A^ , where 

A:=(A^\A^,A^):=a i (G)+a(G), a*(G) := F l c « <8> a S (G x ) d 3 x . 

The physical choice of the coupling function G x = (C?x , G x 2 \ G^) is 

GJk) := -e 1{|k| ll } e -' k - x £ (fc), (2.2) 

for x € R 3 and almost every fc = (k, A) e IR 3 x Z 2 . The parameter A > 
is an ultraviolet cut-off and e G R. (In nature e 2 ~ 1/137 is Sommerfeld's 
fine-structure constant which equals the square of the elementary charge in our 
unitf0-) The values of e and A can be chosen arbitrarily in the whole article. 
Writing 

k x := (fc< 2 ) , -fcW , 0) , k=(fc« fc( 2 \fc( 3 ))eR 3 , (2.3) 
the polarization vectors are given as 

e(k, 0) = jj^r , e(k, 1) = p A e(k, 0) , (2.4) 

for k € R 3 \ {0} with k^ 7^ 0. It is sufficient to determine e almost everywhere. 
Many of our results and estimates do not depend on the special choice of G x . 
If we consider a larger class of coupling functions at a certain point in this 
review we shall explain the required properties of G x at the beginning of the 
corresponding (sub)section. 

2.3 The semi-relativistic Pauli-Fierz and no-pair Hamil- 
tonians 

In order to define the no-pair and semi-relativistic Pauli-Fierz operators we 
recall that the free Dirac operator minimally coupled to A is given as 

3 

D A := a • (-zV + A) + /3 := £ a 3 ( - id Xj + a f (G«) + a(G&)) + f3 . (2.5) 

A straightforward application of Nelson's commutator theorem shows that Da 
is essentially self-adjoint on see, e.g., [32] [40]. We denote its closure starting 
from @> again by the same symbol. As a consequence of (|2.1|) we further have 

D 2 A = Ta © Ta , Ta := (t • HV + A)) 2 + 1 , (2.6) 

1 Energies are measured in units of mc 2 , m denoting the rest mass of an electron and c the 
speed of light. Length, i.e. x, are measured in units of h/(mc), which is the Compton wave 
length divided by 2tt. fi is Planck's constant divided by 2ir. The photon wave vectors k are 
measured in units of 2ir times the inverse Compton wavelength, mc/ti. 
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on 9. In particular, D\ 1 on 9, and since D A is essentially self-adjoint on 
9 we see that ||(D A - z) V>|| ^ C 1 ~ N) Ml> ^ € T>(D A ), zeC, whence 

ct(D a ) C (-oo,-1]U[1,oo). 

Contrary to the usual convention used also in the introduction we define the 
semi-relativistic Pauli-Fierz operator as an operator acting in Jif a priori by 

H^p := (\DA\--fc+H t )<p, <f£3>, (2.7) 
H v F cp := (\D A \+V + H{)tp, y£@. (2.8) 

We shall impose appropriate conditions on the general potential V £ Lf oc (R 3 , R) 
later on. In fact, the operator defined in (|2.7p is a two- fold copy of the one given 
in (II. 4[) since |Da| = Pa^ 2 © 7^ 2 by (|2.6|) . We prefer to consider the operator 
defined by (|2.7p to have a unified notation throughout the following sections. 
The no-pair operator acts in a projected Hilbert space, , given by 

Jlf+ := P+Jf, P+ := l^) Pa), Pa := 1 " ■ (2-9) 

A priori we define it on the dense domain Pj£ £P C > 

i?7 ^+ := P+ (£» A - ft + fff) ^ + eP A + », (2.10) 

F^ p V + := P^(Z>A + F + F f )^ + , ^+eP+^, (2.11) 

where V € Zf^lR 3 , R) satisfies P(P ) C X>(y). In the above definitions we 
have to take care that the right hand sides are actually well-defined as it is, for 
instance, not obvious that ft P^ £? C M . It follows, however, from the proof 

of Lemma 3.4(h) in [37 that P^ maps 9 into V(D ) n T>(Hjf), for every v > 0, 
so that the definitions (|2.10[) and (|2.11j) make sense. 

As soon as we have shown in Section [3] that the quadratic forms of H V F and 
Hy' are semi-bounded below on the dense domains 9 and P^ S>, respectively, 
we may extend them to self-adjoint operators by means of a Friedrichs extension. 
As already mentioned in the introduction there will, however, be critical values 
for 7 above which the quadratic forms are unbounded below in the case of the 
Coulomb potential. We prove in Section [3] that, for H FF , this critical value is 
not less than 

7 C PF := 2/tt, 

which is the critical constant in Kato's inequality, (2/7r)|x| _1 ^ \/—A. In the 
case of the no-pair operator we prove the semi-boundedness of the quadratic 
form of P" p , for all 7 £ [0,7" p ), where 

7c np := 2/(2/^ + ^/2). (2.12) 

The instability of both models above the respective critical values for 7 is shown 
in [281 129) by means of suitable test functions that drive the energy to minus 
infinity. For the definition of iP! p in the case 7 = 7" p see [29] . 
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It has been shown in (T3] that the quadratic form associated with the (elec- 
tronic) Brown- Ravenhall operator, 

:= P+ (Do - ^ P+ , (2.13) 

is semi-bounded on P+ C$°(R 3 , C 4 ) if and only if 7 7° p . Thus, for 7 s$ 7^, it 
has a self-adjoint Friedrichs extension which we again denote by B^ 1 and which 
actually satisfies Tix' inequality, B^ > (1 - 7)P + > 7 € [0,7 c np ]; see gS]. We 
exploit the semi-boundedness of B^ 1 in the proof of Theorem 13.61 below. 

Finally, we introduce a convention used throughout this review: We will 
frequently use the symbols H v , Hfy, 7I, etc. when we treat both the relativistic 
Pauli-Fierz and no-pair operators at the same time; that is, jj is PF or np. 

2.4 How to deal with the non-local terms 

Although general strategies to prove the existence of ground states have been 
developed in the framework of non- relativistic QED [U [6j [17] the application of 
these ideas to the models discussed in this review poses a variety of new mathe- 
matical problems. This is mainly due to the non-locality of the operators |-Da| 
and Pj£ appearing in Hy. In this respect the the no- pair operator is consid- 
erably more difficult to analyze than the semi-relativistic Pauli-Fierz operator 
since also the projected potential and radiation field energies become non-local. 
As a consequence a variety of commutator estimates involving \D\\, Hf, 
cut-off functions etc. is required for a spectral analysis of Hy. Most of these 
commutator estimates are based on the observations and facts we collect in this 
subsection. We shall only present one proof in the present subsection in order 
to illustrate some simple ideas and defer other technical arguments to Section 
We introduce a general hypothesis on the coupling function which is sometimes 
used in the sequel: 

Hypothesis 2.1 The map R 3 x (R 3 x Z2) 3 (x, k) n- G x (fc) is measurable such 
that x 1 — ^ G x (fc) is continuously differentiate, for almost every k, and 

G x (-k,A) = G x (k,A), x e R 3 , a.e. k , A e Z 2 . (2.14) 

There exist <i_ 1 , do , d\ , di € (0, 00) satisfying 

2ju(kr\\G(k)\\ldk^dj, 2 | Jl^gpi^ dfc < df, (2.15) 

where HG^Uco := sup x |G x (fc)|, etc. 

We remark that, if (|2~T4)) is fulfilled, then [A& (x), (y)] = 0, for all j,k € 

{1,2,3}, x, y E R 3 . For later reference we also recall the following well-known 

1/2 

relative bounds, valid for every tp e T>(H { ' ), 

||a-a(G)?M| 2 < d 2 _ x \\Hf /2 *jj\\ 2 , (2.16) 
II" • at(G) V>|| 2 < d\ WH^f + d 2 U\\ 2 - (2.17) 
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In order to cope with the non-locality of P A we write 

Ra(iv) ■= (DA-iy)- 1 , y G R, 

and use the following representation of the sign function of Da as a strongly 
convergent principal value (see Lemma VI. 5. 6 in [23 ), 



5 A <p := -Da I^aI" 1 V = I™ / #a(^)¥>— , V € X. (2. 

T— >oo J TT 



18) 



In addition we observe that 



|£>a| = S A D A , p+ = ll + ls A . (2.19) 

These formulas reduce computations involving \D A \ or P£ to computations in- 
volving Da and integrals over its resolvent. To study the exponential localiza- 
tion it is hence useful to recall that, for all y € R, a € [0, 1), and F € C°°(R 3 , R) 
having a fixed sign and satisfying |VF| ^ a, we have «y € g{D A + ia ■ VF), 

: = e F R A (iy)e- F = {D A + ia-VF-iy)- 1 \ v(e - n , (2.20) 

and 

Pff(<V)ll < - 7 ^=- r ^. (2.21) 
a/I + 1-a 

For classical vector potentials this essentially follows from a computation we 
learned from [9]; see also [38] where (|2.20[) and (|2.21[) are proved in the form 
stated above. It is, however, clear that the arguments in |38] work for a quan- 
tized vector potential, too. Moreover, it is easy to verify that 

[R A (iy), X e F ]e- F = R A (iy) ia ■ (V X + X VF) , (2.22) 

where \ € C°°(Rx; [0, 1]) is some smooth function of the electron coordinates 
and F is as above. Finally, we note that 

||ia-(Vx + xVF)|| < HVxIloo+o, (2.23) 

since \\a ■ v|| = |v|, v 6 R 3 , by the Clifford algebra relations, and |VF| ^ a. 
As an example we treat some commutator estimates whose proofs make use of 
these remarks and a few further useful observations. 

Lemma 2.2 Assume that G x fulfills Hypothesis \ 2. 11 Let x and F be as above, 
assume additionally that F is bounded, and set Hf := Hf + E, for some suf- 
ficiently large E ^ 1 (depending on d\). Let V G L 1 1 oc (R 3 ,R) be relatively 
form-bounded with respect to y/— A. Then, for all ao,K G [0,1), v ^ 0, and 
a G [0,oo], 

|| |D A | K [F+, xe F ]e- F || < const(a ,K)-(a+||Vx||oo), (2.24) 
H^n-PA.X^e^^ni < const(oo,i/)-(o+||Vx||oo), (2-25) 
||l^r /2 [F+, X e F ] e - F iJ- 1/2 || < const( ao ,^)-(a + ||Vx||oo). (2.26) 
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Notice that the ao-dependence of the constants originates from the singularity 
at a = 1 in (|2 . 2 1 1) . Notice also that we may choose V — |x| _1 in (|2.26p in view 
of Kato's inequality. 

Before the proof we further remark that all operators appearing in the norms 
in (|2.24p - (|2.26p and in similar estimates below are always well-defined a priori 
on & and have unique bounded extensions to the whole Hilbert space. In fact, 
Pj[ 9 C V(D ) n n„> v ( H i) as we have recalled from [37] above already. To 
simplify the presentation we shall not comment on this anymore from now on. 

Proof: We use the fact that an operator, T, acting in some Hilbert space 
is bounded if and only if sup^n \\^p\\ = i \(<p \ Ti/j)\ is bounded in which case it is 
equal to \\T\\. Here it is sufficient to take the supremum over all normalized ip 
and ip from a dense set which is a core for T. Combining (|2.18[) . (|2.19p . and 
(I2.22p we find, for all normalized tp,ij) € 

(\D A \^\[P+, X e F ]e- F ^) 

= lim / (\D A \ K ip\R A (iy)ioc-(V X + X^F)R F (iy)^)^-. 

T^°°J~T 27T 

On account of || \D A \ K R A (iy)\\ ^ const(/s)(l + j, 2 )-i/2+ K /2 ; ([^1) and (|2~231 
we see that the scalar product under the integral sign defines some Lebesgue 
integrable function of y and 



{\D A \ K <p\[P£, X e F ]e- F il))\ <const(K)(||Vx||oc+a) 



dy 



( 1 + y 2)l-«/2 



where the last integral is finite. Therefore, [Pj^ , x eF ] e F ip belongs to the 
domain of (\D A \ K )* = \D A \ K and the first bound ([234]) follows. 

In order to prove the second bound (|2 . 25|) we introduce another little tool 
which turns out to be useful in our whole analysis. Namely, if E 1 is suf- 
ficiently large depending on di we can construct T F (iy) £ Jif(Jf?) such that 
R A (iy) H7 V — H7 V R^(iy) T F (iy), for every y £ R, and such that the norm of 
T F (iy) is uniformly bounded with respect to y £ R; see Corollary 19.21 below. 
(In particular, R^(iy) maps T>(H") into itself.) Therefore, 

\(H?<p\ R A (iy) ia ■ (V X + X ^F) R F (iy) fif $ ) \ 

= \(<p\R A (iy) T° (iy) ia-(V X + xVF)R F A (iy) T F (iy) j> ) \ 
< C(\\V X \\oo+a)(l + y 2 r\ 

and it is clear from the argument above how to derive (|2.25p . 

The last bound (12.26)) follows from the first two and the inequality |V|/C7 ^ 
I -Da | + H{ + E proved later on in Theorem EOl □ 



3 Self-adjointness 

As it is obvious from the definitions in the preceding section the operators Hq 
are positive. In this section we present some basic relative bounds that allow to 
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define the perturbed operators H v as self-adjoint Friedrichs extensions. As a 
rule we denote the self-adjoint extensions of H v or H$ - which are only defined 
on St and P A 2$ so far - again by the same symbols. For suitable V, we are 
also able to characterize the quadratic form domains of H v which turn out to 
be the spaces of all vectors with finite kinetic and radiation field energy; see 
Theorems 13.41 and 13.61 below. 

Before we present the afore-mentioned results we discuss various (essentially 
well-known) diamagnetic inequalities in QED; see Theorem 13.21 below. Since 
these estimates are of independent interest we decided to present one way to 
derive them (adapted from [47] ) in detail which has not been worked out in the 
literature before, as it seems to us. 

3.1 Diamagnetic inequalities in QED 

In this subsection it is sufficient to assume that 



where R 3 x (R 3 x Z 2 ) 3 (x, k) H> G x (fc) is measurable such that x h-> G x (fc) is 
continuously differentiable, for almost every k, and 



The following result is probably well-known but the argument sketched in its 
proof might be new. 

Lemma 3.1 Let A 0. Under the above condition on G x the operator (— iV + 
A) 2 + \H{ is essentially self-adjoint on 2>. 

Proof: It is a standard exercise to show that {— iV, A} + A 2 is a small 
operator perturbation of — A + cHf, provided that c > is chosen sufficiently 
large depending on G x - In particular, N := (— iV + A) 2 + c Hi + 1 is essentially 
self-adjoint on any core of — A + cH{ and, in particular, on Si . In the next step 
we apply Nelson's commutator theorem with the closure of M starting from S 
as test operator to conclude. □ 

We denote the closure of (— iV + A) 2 starting from by ta- For every 0, ip 6 
= L 2 (R 3 ,C 4 ® J^b), we write (<f>\ip) for the (partial) scalar product on 
C 4 J^b and denote [y](x) := (y(x) | ^(x) ) 1 / 2 . Furthermore, we set 5^(x) := 
0(x), for 0(x) ^ 0, and ^(x) = 0, for 0(x) = 0. 






Theorem 3.2 (i) Let (f> G V(t a ). Then {tg G ff^lR 3 ), and 

(ry | -AW) lS(R3) < Re f ^(S^xJKta^Kx)))^ 



J K 3 



(3.1) 
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for all rj G iJ 1 (R 3 ) 7 77 0. in particular, for r\ = [(/>], 

(Ml -Ai) iS(R3) < ^|r A 0). (3.2) 
fuj Let <p G V{t^ 2 ). Then [0] G i? 1 /2( R 3^ and 

(tjl-Z^AM^CRS) < Re f V (x)(S^x)\(r^ 2 ^(x))d 3 x, (3.3) 

JR 3 

for all rj G i/ 1 / 2 (IR 3 ) 7 77 0. In particular, for rj = ftfij, 

([0]|V=AM)l> (R .) < (0|rI /2 0). (3.4) 

fraj For aZ/ V G an ^ ^ G [0>°°)> we have, almost everywhere on R 3 , 

[e-'^J ^ e-*(- A )[^], (3.5) 

\e~ tT - 2 n sc e-^m- (3-6) 

Remark 3.3 (1) Arguing as in Theorem 7.21 of [30] with the corresponding 
changes as in the proof below one can easily verify that, for <j) £ @, we have 
[01 G H\R 3 ) and |V [0](x)| [(-iV + A)0](x), for a.e. x G R 3 . 

(2) In [19) diamagnetic inequalities for infra-red regularized vector potentials 
have been proved by means of dressing transformations. For an alternative 
proof using functional integrals see [20] . If all components of the vector poten- 
tial commute, [A^ (x), (y)] = 0, then one can also reduce the diamagnetic 
inequalities to classical ones by diagonalizing all components A^'(x) simulta- 
neously; this argument due to J. Frohlich is mentioned in pQ. The proofs given 
here are variants of the ones presented in [43] [47] . 

Proof: Let e > 0. First, we assume that <j> G @ and set u e := \/PP + e 2 G 
C°°(IR 3 ,[R). Since A^(x) is symmetric on for every x G R 3 , we have 
Re ( <p(x) I i A {j) (x) (7>(x) ) = 0, thus 

u E Wu E = - V u\ = Rc(4>\W<p) = Re ( <f> | (V + iA) 4> ) , (3.7) 
In particular, 

|V u e \ < M [(V + iA) 0] sc [(V + iA) 0] . (3.8) 

u e 

Taking the divergence of Q3.7P we obtain 

I V u E \' 2 + u £ A u £ = Re ( | (V + iA) ) + Re ( | V (V + iA) ) 

= [(V + zA)0J 2 -Rc (0|r A 0). (3-9) 
Combining this identity with (13.81) we arrive at 

- Au e < Re (uJVIta^). (3.10) 
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Now, assume that </> £ V{ta)- Since ta is essentially self-adjoint on 3) we find 
<p n £ n £ N, such that — > <p and T\<p n — > ta</> in c^ 5 . On account of 
PTTU]) we have 

(-A7 7 )(x)4 n )(x)d 3 x < Re<r ? ( U ("))- 1 0„|rA^}, (3.11) 

R 3 

for all Schwartz functions r? £ .Y{R 3 ), r) > 0, where 4"' := v^PnP +^ 2 : n £ M. 
Passing to appropriate subsequences if necessary we may assume that [0„J — > 
[01 and, hence, t4"' ) — » u 6 almost everywhere. Using that u^ 1 , (■t4"' ) )~ :l $C 1/e, 
it is easy to see that rj (ite™^) -1 0„ — > 7/7i~ 1 in By virtue of the Riesz- 
Fischer theorem we further find a square-integrable majorant for the sequence 
([0„J). We can thus pass to the limit n — > oo in p. Ill) to get, for all 77 £ Y(R 3 ), 
r) ^ 0, and <j) £ 2?(ta), 

/ (-Ar;)(x) U£ (x)d 3 xs;Re / r?(x) ( u" 1 (x)0(x) | (r A 0)(x) ) d 3 x. (3.12) 

JR 3 JR 3 

Here we may take the limit e — > by means of the dominated convergence 
theorem (with the majorant 77 [ta0] on the right hand side) to obtain, for all 
77 £ y(R 3 ), 77 ^ 0, and <p £ V(t a ), 



(-A77)(x)[0Kx)d 3 x sc: Re / r?(x) ( S*(x) | (r A 0)(x) ) d 3 x . (3.13) 

R 3 JR 3 

Adding /»7A[0J with A > to both sides we obtain 

[(-A + A)77](x)[01(x)rf 3 x^ Re / tj(x) ( | (r A + A)0)(x) d 3 x 



R 3 

^ / 7,(x)[(TA+A)0J(x)d 3 X. 

JR 3 

Let < x e C^ (R 3 ) and V € ^ Since (-A + A) -x , A > 0, is positivity 
preserving, we may then choose 77 := (—A + A) _1 x £ =5^(R 3 ) and <p '■= ( t a + 
A) _1 7/> £ 2?(ta) and arrive at 



X (x)[(r A + A)-V](x)d 3 x < / x( x ) [ (—A + A) [ V>] ](x) crx . 

R 3 JR 3 

Since \ € Qf([R 3 ) is arbitrary we find [ (ta + A) _1 i/> J ^ (-A+A) _1 [^], almost 
everywhere on R 3 , and by induction (see 47 , for the same argument) we get, 
for all 77 £ N and £ > 0, 

[(tr(TA+t)- n v>] ^ (tn-A+trw 

Both sides converge almost everywhere along some subsequence to [e - * TA 7/>] and 
e~*( _A )[7/)J respectively, and (|3.5[) follows. Equation (|3 . 6[) follows from ()3.5|) . the 
spectral calculus, the subordination identity 

e^ 1/2 = fe— *"V(4.) * , t.A^o, 
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and the properties of the Bochner-Lebesgue integral. In the remaining part of 
the proof we derive (following again [47]) (|3 . 1 1) and (|3.3I) at the same time. To 
this end let v G {1/2, 1} and <f> G £>(t a ). 

On account of ( cj) | e _ * T A^) < / [0] [e"* T A0] Equations J33]) and J3^J imply 

(^le-^A^X fM e - ( (- A )"W, thus 



(^It-Hl-e-^) > / r^l-e^lMf (e)d 3 e- 

Since G 2?(t"a) the hmit i — > exists on the left hand side of the previous 
inequality. By the monotone convergence theorem we conclude that the limit 
t — > of the right hand side exists, too, and 

Jr 3 

Hence [0J £ H"(R 3 ) and (j3~2"]l and ([3~3]) hold true. Using this we may take the 
derivatives at t = on the left and right sides of the following consequence of 
fiSfl and (BSJ, 



Re / 7 ? (x)(5 |e- tr A0)(x)d 3 x^ / 7?(x)[ e - tT A 0](x) d 3 x 

y R3 y R a (314) 



to get -Re (t7S |t£0) < -( V \ (-A) v [^]) £a(RS) , for all r? G i^(R 3 ),r/ ^ 0. 
Here we have also used that all expressions in (I3.14[) are equal to /^[^] at 
* = 0. □ 



3.2 Semi-boundedness 

The following theorem is a slight generalization of a result from [37]. Its proof 

is based on two basic steps: The first one follows immediately from the diamag- 

netic inequalities by means of which the form bounds the potential satisfies by 

assumption can be turned into form bounds with respect to the scalar operators 
l li 

t a or t a- After that we use relative bounds on the magnetic field to include 
spin, t a ^ D\ + + const. To complete the square on the right hand side 
of the previous bound we employ the inequality p.20[) below. After completing 

1/2 

the square we take square roots on both sides to obtain a bound on t a . 

Theorem 3.4 For v G {1/2, 1}, let V v G L\ oc {¥?, R) and assume that there is 
some c v such that, for every ip G H U (R 3 ), 

(<p\V„(p) L 2 {R3) < (tp\ (~A) U <p) L 2 (R 3 } +c„ || <f\\h(R3)- (3.15) 

Then there is some C G (0,oo) such that, for all G x fulfilling Hypothesis \2. 11 
4> G @, and S > 0, 

{<t>\V v <t>) < (4>\ (\D A \ + 6H l + (6- 1 + C6)dfj 2u cp) + c u Uf. (3.16) 



17 



In particular, 

i|||xrV|| 2 < || (\D A \+SHf + (S- 1 +CS)dj)cj ) \\ 2 , (3.17) 
for all and 

-A < \D A \ + SH [ + (S- 1 + CS)d 2 1 , (3.18) 

7T |x| 

in the sense of quadratic forms on '3. Therefore, Hy F and H PF , 7 £ [0,7,? F ], 
have self-adjoint Friedrichs extensions - henceforth again denoted by the same 
symbols - and 3 is a form core for these extensions. Moreover, for a £ [0, 1) and 
7 G [0, 7c PF ), we know that Q(H™ /a ) = S(i?? F ) = 2(# PF ) = Q(\D \)nQ(H f ). 

Proof: First, we show that Q(H PF ) = Q(\D \) n Q{H f ); see (SUES]. The 
remaining statements on form domains will then be a consequence of (|3.16j) and 
(|3.18p . In fact, this follows from the bounds [37] 

HiA)r /2 (SA-soW 1/2 MC; iiic A r /2 (^A-5o)i?r 1/2 Nc, 

where Hf := iff + E with E ^ 1 + (4di) 2 . (These bounds are derived exactly as 
in Lemma IOI below.) Together with |D A | - |D | = D (S A - S ) + a ■ AS A 
and (|^T6]l fc (P371) the first bound implies 

|< <p I (|D A | - \D \)<p)\ < C || \D \V 2 <p\\ H^Vll 

<C"(^|(|D |+fff)^}, (3.19) 

for all Lp £ 3>. Analogously, the second bound implies (13.19[) with \Dq\ replaced 
by \D A \ on the right hand side. Consequently, the form norms of |D A | + iff 
and I D 1 + Hi are equivalent on <® which implies Q(H PF ) = Q(|D |) n Q(H { ). 

All details missing in the proof of (|3.16|) - (|3.18p sketched below can be found 
in [37] . We set Hf := iff + E, for E > 0. Besides some standard arguments the 
main ingredient in this proof is the following bound proven in [371 Lemma 4.1]: 
We find some constant, C > 0, such that, for all E > C d\ and <j> 6 S>, 

Re(|D A |</>|fff </>) > 0. (3.20) 

This estimate follows from the following identity Re (|D A | Hf) = H { ' (|D A | - 
T)H\ 12 on 9, where 

T := Re{[|D A |, Hf 1/2 ]Hl /2 } ^ e |D A | + e" 1 const di/f±: 1/2 , 

for £ G (0, 1] and -E (4c?i) 2 , as we shall see at the end of Subsection 19.21 To 
make use of the bound (|3.20[) we recall that, since [A^ (x) , A^ k ' (y)] = 0, we 
have 

D A <t> = ta4> + + 4>, 4>£3>, (3.21) 
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where the entries of the formal three-vector S are Sj — <7j <X> I2 and B is the 
magnetic field, i.e. S • B = S • at(V x A G) + S • a(V x AG). By (j2~T6]) with 
(V x A G, di) instead of (G, we have, for all S > and <f) € @, 

|< I S • B 0>| sC 2di||^|| ||# f 1/2 </>|| < 8{4>\{H S + d- 2 d\) 4>) . (3.22) 

Choosing E = (S- 2 + C) d\ we infer from (l3~2Ul) - (l3~22l . for all (j> e S>, 

(4>\ta^) < < D A <j>\ D A <j>) +5(4>\H {( j>) -II0H 2 (3.23) 
< (Dj^cj)] D A <j)) + (<b\5 2 H 2 0) + 2Rc (\D A \<t>\5H {( j>) 

= ||(|£> A |+<Sff f )0|| 2 . 

Furthermore, since the square root is operator monotone it follows from (|3 . 23[) 
that ( <p J Tj[ 2 </>) < ( </> I (I -Da I + 5 H{) cf>). Using the diamagnetic inequalities 
(j3~2l) and p~3} we further find, for v € {1/2, 1}, 

= ([^]|K[0])L 2 ( R 3) ^ ([0]|(-A)-[0]) L2(R 3 )+Cl ,||0|| 2 

< (<\>\rt 4>)+c„U\\\ (3.24) 

and we conclude that Q3.16P holds true. Inequalities (|3.17p and (|3.18p follow 
from p. 241) together with Hardy's and Kato's inequality, respectively. □ 

Theorem 13.41 has a straightforward extension to the case of N electrons [35] . 
We discuss this extension in the next corollary mainly since its proof gives the 
opportunity to introduce some identities and estimates which are used later on. 
Let and 3>n, N £ N, be defined in the same way as Jf? and *2) but with the 
I 2 (R 3 ,C 4 )=i 2 (K 3 xZ 4 ) replaced by L 2 ((R 3 x Z±) N ). The spatial coordinates 
of the i-th electron are denoted by x< £ R 3 and we designate an operator acting 
only on Xj, the i-th spinor components, and on the photon field by a superscript 
(z). 

Corollary 3.5 Assume that N,K e N, e > 0, 71, •••,7k G (0,2/7r], and 
{Ri,...,R K } C R 3 . Then 

N N K 2 

i=l i=l fc=l 1 1 i<j 1 J 1 

/or every S > 0, in £/ie sense of quadratic forms on &n ■ 

Proof: In view of (I3.18|) we only have to explain how to localize the non- 
local kinetic energy terms. To begin with we note the following special cases of 
(|2~24l and (|9T0]) . respectively: For every X € C°°(R X , [0, 1]), 

II [x, S A ] || s$ || Vxlloo , \\D A [ X , [X, 5 A ] ] II < 2 HVxIlL • (3-26) 

Let B r (z) denote the open ball in R 3 of radius r > centered at z G R 3 . 
We set g := min{|Rfe — R^| : k ^ £}/2 and pick a smooth partition of unity 
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on R 3 , {xk}k =Q , such that Xk = 1 on B g / 2 (R k ) and supp(xk) C B e (Rfc), for 
fc = 1, . . . , K, and such that ^fc=o Xfe = 1- Combining the following IMS type 
localization formula, 

K 1 
|£>a| = J2{xk\DA\Xk + 2 [xfc, lXk,\DA\]] } on#, (3.27) 

fc=0 

and the identities 

[Xk,\Xk,\D A \]] = 2ia-(Vxk)[Xk,SA]+D A [ Xk , \xk,S A ]] (3.28) 
and ||a • Vxfc|| = |Vx/c| with (|3.26[) , we obtain 

|| [Xk, [Xk,\D A \]] || < 4||Vxfc||L. ( 3 - 29 ) 

for all k G {0, . . . , K}. Since we are able to localize the kinetic energy terms 
and since, by the choice of the partition of unity, the functions R 3 9 x n- 
|x- Rfc| _1 xK x ) are bounded, for k G {1, . . . , K}, I G {0, . . . , K}, k ^ £, the 
bound (|3.25l) is now an immediate consequence of (|3.18[) . □ 

Next, we discuss the semi-boundedness of the no-pair operator. The idea is to 
reduce the stability of the no-pair operator to the one of the purely electronic 
Brown- Ravenhall operator. 

Theorem 3.6 Assume that G x fulfills Hypothesis \2.1\ Let 6 G (0, 1] and let 

V G Lj 2 oc (lR 3 , R) be form bounded with respect to \J —A and satisfy 

(P+ l p\VP+^) ^a( lf \D P+^)+b\\P+ l p\\ 2 , V ECZ°(R 3 ), (3.30) 

for some a G (0,1) and b ^ 0. Then there exist constants cy,C G (0, oo), 
C = C(S,V,d-i,di), such that, for all (p+ G P A @, \\<p + \\ = 1, 

(<P + \(D A + V + SH t ) <p + > cy( tp+ | \D \ <p + )-C, (3.31) 

and in particular, for every 7 G [0, 7" p ), 

(<p+\{D A - ft+6H t )<p + ) > c^(<p + \\D \<p + )-C(5,~/,d-i,d 1 ). (3.32) 

Therefore, H 1 ^ and H^ p have self-adjoint Friedrichs extensions - henceforth 
again denoted by the same symbols - and P^ & is a form core for these exten- 
sions. Furthermore, Q(H* P ) = Q(H^) = Q(Hq P ) = Q(|A>|)nQ(fff )nRanP+. 

Proof: The statement on the form domains follows from Theorem 3.2 of 
[29] . (See also Section 3.4 of the first preprint version of [29] available on the 
arXiv for an alternative proof.) The estimate (|3.32[) is derived in [37] and we 
shall outline its proof in what follows. 
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We pick some p > 1 with p a < 1 and write, for ip + £ Pj£ fF, 

( ^+ | (D A + V) V +) = p-\ V + | P+ (D + P V) P+ ^ ) 
+ (l-p- 1 )(^ + \P+D oV +) 
+ (lp+\ol- Alp+) (3.33) 
+ ( <p+ | P- (Do + 7) P " ^+ ) 
+ 2Re(^+|P +FP -^+). 

The estimate (j3.32|) is based on the identity above and the following bound 
on the difference between the spectral projections with and without field (see 
Lemma [6.31 for similar statements): For E 1 + (4<ii) 2 , there is a constant, 
C = C(d-i,do) > 0, such that 

|||I>o| 3/ *(J^-P|)flr 1/a || < C, (3.34) 

where H f = H { + E. On account of P Lp+ = {P - P A ) <p+, for ip+ £ P£ @, 
and (|3.34p we have, for every e > 0, 

|| |D | 1/2 Po ^ f < II I A,| 1/4 P ^+ || || |Po| 3/4 (P Pa) V + II 

^ 1 II I n 1 1/2 d- + II 2 i C(£,d_i,do) || D - + i|2 £ || fr l/2 + 1 1 22 

< 2 II l-°o| 7 P ¥> || + ^ || p o f || +2 II ff f P II ' 

that is, 

HIA)! 172 ^"^!! 2 < e\\Hy 2 p + \\ 2 + C(e,d„ 1 ,d )\\P p+ \\ 2 . (3.35) 
By virtue of |V| ^ C \Do\, the previous estimate further implies, for every r > 0, 
|(P +^+|^P -^+)| 

< r || |D | 1/2 P + f + £ || P f 1/2 ^+ f + C £ , T || P - ^+ f. (3.36) 

Here the second term on the RHS of (|3.33|) can be used to control the first 
term on the RHS of (|3.36p . Recalling the definition (I2.13[) and applying (|2.16p . 
(|3.35p . (|3.36p . and (|3.30p to the various terms in (|3.33p we thus find, for every 
S £ (0,1], some constant, C = C(6, p,d-t,di) £ (0, oo), such that 

(<P + \(D A + V + 6H { ) <p+ ) > c a ^+ | D P+ <p+ ) - C || tp + 1| 2 . 

Using (|3.35p once more to replace Do Pq by \Dq\ on the right hand side, we 
arrive at the first asserted estimate. According to the remarks made below 
(|2.13p the first estimate applies in particular to the Coulomb potential, as long 
as 7 e [0,7° p ). □ 
From the previous theorem and our commutator estimates one can also infer the 
semi-boundedness of a no-pair operator for N £ N electrons and K £ N static 
nuclei, analogously to Corollary 13.51 as long as all Coulomb coupling constants 
7i, ... ,7k are less than 7" p ; see Proposition A. 2 of [3"6] . 

Since we are addressing the question of finding distinguished self-adjoint 
realizations of PjJ it is also natural to state the following theorem whose proof 
can be found in Corollary 3.4 of [2"5] . 
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Theorem 3.7 Let 7 e [0,1/2) and assume that G x fulfills Hypothesis lff.il 
XTien -ff PF and i?" p are essentially self-adjoint on and P£ 3, respectively. 

For sufficiently small values of |e| and/or A, the essential self-adjointness of H^ F 
has been shown earlier in [40] . 

4 Bounds on the ionization energy 

As a first step towards the proof of the existence of ground states we need to 
show that binding occurs in the atomic system defined by Hy in the sense that 
inf cr[-f/y] < inf cr[-ffg]. This information will be exploited mathematically when 
we apply a bound on the spatial localization of low-lying spectral subspaces of 
Hy from [37]. The localization estimate in turn enters into the proof of the 
existence of ground states at various places, for instance, into the derivation 
of the infra-red estimates and into the compactness argument given in Subsec- 
tion 18.11 Theorem 14.11 below is the main result of this section. In its statement 
we abbreviate (ft € {PF,np}) 

E v := inf a[H v ] , E\ := inf a[H\] , 7 e (0, 7 ») , E» := inf a[H*] , 

where V satisfies the conditions under which H v has been defined in the pre- 
vious section. To simplify the exposition we only consider the physical choice 
of the coupling function G x given in (|2.2I) . as always for arbitrary values of e 
and A. Our proofs work, however, equally well for other coupling functions, 
for instance, for the infra-red cut-off and discretized coupling functions intro- 
duced in Section [51 and we obtain uniform bounds on the binding energies in 
these cases. If we consider coupling functions other than (|2.2p then the unitary 
transformation U introduced below has to be changed accordingly; see [351 US] • 

Theorem 4.1 (i) Let V € L^ oc (lR 3 ,lR) be form bounded with respect to \J — A 
with form bound less than or equal to one. (So V fulfills (|3.15|) with v = 1/2.) 
Define the self-adjoint operator hy '■= — A + V by means of a Friedrichs 
extension starting from C^°(IR 3 ) and assume that vaia[hv] is an eigenvalue. 
Then 

S pf„£PF ^ l-M a [hy]. (4.1) 
In particular, for 7 € [0, 7 PF ], 

S PF - £ PF ^ 1 - inf a[VT^A - ft] . (4.2) 

(ii) Let V be relatively form bounded with respect to \J — A and assume that 
V satisfies (|3.30[) with a e (0, 1). Additionally, assume there exist r ^ 1, c > 0, 
and 8 £ (0, 2) such that 

y(x)<-c|x|^ 2 , |x|>r. 

Then 

E np _ E n P > q ^ ^3) 
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and in particular, for 7 £ (0,7,? p ), 

Remark 4.2 (i) The bound (|4. 1[) has been obtained first in 21 (under the 
assumption that H FF be essentially self- adjoint which, in the case V = — A, is 
true, at least for all 7 < 1/2). The result of [21] improved a lower bound on 
the binding energy in an earlier preprint version of |28j . The latter was given 
in terms of the non-relativistic ground state energy of an electronic Schrodinger 
operator. 

(ii) In a forthcoming work of the first two authors [27] is it shown that the 
inequalities (I4.1j) and (|4.2[) are actually strict, for all e, A > 0. Moreover, there 
is a certain class of short-range potentials V such that E PF — E FF > and 
in particular - according to the present article - E V F is an eigenvalue of H V F 
although infer [hy] — 1 and 1 is not an eigenvalue of hy. This effect is called 
enhanced binding due to the quantized radiation field and we are again able to 
prove its occurrence, for arbitrary large values of e and A. There are numerous 
results on enhanced binding in non-relativistic QED; up to now complete proofs 
were, however, available only for small e. The proofs in |27j extend the ideas 
and methods underlying the proof of Theorem 14.11 given below. 

Because of lack of space we shall only describe the proof of Theorem 14.11 
for the semi-relativistic Pauli-Fierz operator [28] in detail. The proof of (|4.3I) 
follows similar lines and can be found in [29]; see also Remark (|4.3[) below. 

Our proof of (|4.ip and (|4.3|) is based on a direct fiber decomposition of ffl 
with respect to fixed values of the total momentum p <g) 1 + 1 ® pf, where 
p := — iV x and 

p f := dT(k) := (dr(fc«), dT(k^) , dT(k^)) (4.4) 

is the photon momentum operator. In fact, a conjugation of the Dirac operator 
with the unitary operator e* Pf ' x - which is simply a multiplication with the 
phase (R 3 )" 9 (ki, . . . , k„) e *(ki+-+k„)-x in each Fock gpace sector 
- yields 

e lpf ' x D A e- lpt - x = a ■ (p - p f + A(0)) + , 

and a further conjugation with the Fourier transform, T : L 2 (R 3 C ) — > L 2 (Rp), 
turns the latter expressions into 

(J-(8)l)e jp '' x D A e- jpf - x (J : "- 1 (g)l) - / D{P)d 3 P. (4.5) 

Here the operators 

5(P) := a-(P-p f +A(0))+ J 9, PeR 3 , 

acting in C 4 <£> #b P^] j are fiber Hamiltonians of the transformed Dirac operator 
in (|4.5|) with respect to the isomorphism 

3tf = L 2 ([Rf>,C 4 )<g>^ b [jr] Si / C 4 (8)^" b [jr]d 3 P. (4.6) 
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(In particular, the transformed Dirac operator in (|4.5p again acts in J^, where 
the variable in the first tensor factor, P, is now interpreted as the total momen- 
tum of the combined electron-photon system.) Accordingly, we have the direct 
integral representation (compare, e.g., Theorem XIII. 85 in [4"4"] ) 

(F ® 1) e ipf ' x PF e- jpt ' x (J 7-1 9)1) = l H^TP) d 3 P , (4.7) 

Jr 3 

where 

if PF (P) := |S(P)|+# f . 

Let e > 0. Then we know that the Lebesgue measure of the set of all P G R 3 
satisfying c[i/ PF (P)] n (S PF —e, S PF +e) 7^ is strictly positive. In particular, 
we find some P* € R 3 and some normalized y>* € <2(i/ PF (P*)) such that 

<^|^o PF (P*)^> WbW <S PF + £ . (4.8) 

We define the unitary transformation 

U := e J (Pf- p *)-x (4 9) 

and observe as above that 

UD A U* = 5 P (P*) := a- (p + P*-p f + A(0))+/3. 

It is sufficient to prove the bound (14. ip for the unitarily equivalent operator 

UH(/U* = |5 P (P*)| + K + ff f . (4.10) 

PROOF: (Theorem \4-l\ The semi-relativistic Pauli-Fierz case.) Let e > 
and Pi be as in the paragraph preceding the statement. We abbreviate t* := 
P^ - p f + A(0) and, for r) > 0, 

Ri(v) ■= (P 2 + (a-t*) 2 + T7 + Mv) ■= ((a-(P + t*)) 2 + »7 + l)"*- 

Since the anti-commutator of a ■ p and a • t+ is equal to 2 p • t+ it holds (a ■ 
(p + t*)) 2 = (a • ti) 2 + 2 p • + p 2 . We deduce that, for any ip e @, 

-ifcfa)? = -R 2 (v) [p 2 + (a ■ t*) 2 + 1 + ??] flifajp 

= R 2 {ri)[2 V -U]R 1 {r t )(p-R 1 {ri)ip.. (4.11) 

We use the following formula, for a self-adjoint operator T > 0, 

- r (i--JL\ {p ^L, <peV(T), (4.12) 
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and the resolvent identity (|4.11[) to obtain, for any tp £ £>, 

<^|(v/(a-(p + t,)) 2 + l - y/p* + (a-U)* + l)<p) 

= 1™ {<p\{Ri(ii)-M*i))<p)y/ri — 
Jo n 

= [°° (R 2 ( V ) V >\[2p-U}R 1 ( V ) 9 )^j^ 
Jo 7r 

= f (<p\Ri(ji) [2p-U]Ri(ri)<p)y/v — 
Jo 7r 

- / (tp\R 1 (r 1 )[2p-U}R 2 (r l )[2p-U}R 1 (j 1 )cp)^^ 
Jo 7r 

< f < ^ j iZxCf?) [2 p ■ t*] i2i(t?) ^ >V^— • (4-13) 
Jo 7r 

In the last step we used the positivity of i?2( 7 7)- We consider now ip := ipi ® ip 2 
where ip x £ C£° (R 3 ,R) and ip 2 € C 4 <g) % with ||^|| = 1, j = 1,2. Writing 
$2(£, »?) := (£ 2 + (a ■ t*) 2 + 1 + if}- 1 ip 2 = $ 2 (-£, V) we nnd that 

= [ $-(®2&v)\t*<S>2&v))\Mt)\ 2 d 3 d = 0, (4.14) 

J[R 3 

due to the fact that ipi is real and, hence, |<pi(£)l = Iv'iX - Furthermore, 

|S p (p,)hVMp+W+T, |5cp*)| = V(o-t*) 2 + i. (4-i5) 

Vp 2 + (a ■ t*) 2 + 1 < ^(a ' t*) 2 + 1 + vV + 1 ~ 1 , (4-16) 

where we used [p 2 , (a • t*) 2 ] = in the second line. Combining (|4.13l) ~ (|4.16p 
we arrive at 

( if | UH^ F U*ip) «S (cp 2 | H PF (P+) <^ 2 > + (<Pi | V?>1> - 1- 

By a limiting argument the previous inequality extends to any real- valued ipi € 
Q(/iy). We choose cp\ to be he normalized, strictly positive eigenfunction of hy 
corresponding to the eigenvalue at the bottom of its spectrum and ip% = ip+. By 
the choice of v?+ in (|4.8[) . where e > is arbitrary, this proves the assertion. □ 

Remark 4.3 As already mentioned the proof of Theorem 14.11 for the no-pair 
operator employs ideas similar to those described above. However, due to the 
more complex structure of the no-pair Hamiltonian the resulting bound (|4.3I) 
is not as satisfactory as the one for Hy F . Again, we have the representation 
tfo P -/^o P (P)rf 3 Pwith 

iC(P) = P(P) (D(P) + H f ) P(P) , P(P) := l [0 ,oo)(5(P)) . 
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Given e > 0, we again find P* € R 3 and some normalized <y9* = P(P) <p+ £ 
Q(ff n P(P*)) such that (¥J*|flS P (P*)v*>c*®*i,[jf] < S"P + £ . Conjugating 
H^ p with E7 defined as in (|4.9I) we obtain 

[Tfl?>tf* = P p (P,)(5 p (P,) + ra f )P p (P,), (4.17) 

where P p (P*) := l^.oo) (-Dp(P*))- As a test function we now choose := 
Pp(P*)(XR ® V*); where Xii is some normalized smooth and non- negative func- 
tion supported in {R ^ |x| ^ 2R} with i? ^ 1. It turns out that = 
1 + 0(1/R 2 ). We then exploit the fact that the negative contribution of the 
potential in (|4.17p to the expectation value of decays as l/i? 2_1? whereas 
all other terms yield a contribution S np + e + 0(1/R 2 ), as R tends to infinity. 
To show all this it is convenient to work on a non-projected Hilbert space by 
adding a suitable "positronic" no-pair operator to H™ p (similarly as in (15 . 1 1) 
below). For then one can again make use of the bounds f|4.13[) — f|4.16[) derived 
in the previous proof; see Section V of 29 . (It is entirely obvious that the 
Coulomb potential can be replaced by a more general one in Section V of [29].) 

5 Exponential localization 

Our next aim is to discuss the exponential localization with respect to the 
electron coordinates of low-lying spectral subspaces of the no-pair and semi- 
relativistic Pauli-Fierz operators. Since the multiplication with some exponen- 
tial weight function, e F , acting on the electron coordinates does not map the 
projected Hilbert space Jf A into itself it is convenient to extend the no-pair 
operator to some continuously invertible operator on the whole Hilbert space 
Jrf? in the discussion below. Therefore, we set 

H v p := H$> + H$»'- , := P A (\D A \ + H f ) P A . (5.1) 

In Lemma [8.31 below we show that E np = inf cr[_ffQ P '~] by constructing some 
anti-linear map r : Jf? — > J$? with Hq P t = t H^ P: ~ . Therefore, 

H^ p Si E np . (5.2) 

To unify the notation we further set Hy F := Hy F and write H v , S", etc., when 
we treat both Hy F and H v p at the same time. In the whole section we assume 
that G x fulfills Hypothesis ED 

Theorem 5.1 (i) Let V £ if oc (IR 3 , R) be relatively form bounded with respect 
to \] — A with relative form bound less than or equal to one. Moreover, assume 
that 

3r>l: sup |V(x)| < oo and V(x) |xK °° > o, (5.3) 

|x|>r/4 

Define 

p(a) := 1 - (1 - a 2 ) 1 / 2 , a £ [0, 1) , 
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and let I C (—00, S PF ) be some compact interval. Then there exists k £ (0,oo), 
such that, for all a £ (0, 1) satisfying e :— E PF — sup/ — p(a) £ (0, 1], we have 

|| e a ^ || < k (E PF - E V F ) e fe / £ . (5.4) 

(ii) Assume that V £ Lf oc (R 3 , R) satisfies i? 1 (K 3 ) C X>(V) (Ww'c/i implies \V\ < 
const IDoIJ as we/Z as (pT3T)l) a < 1 and ([53)1 . iei / C (-00, £ np ) be 

some compact interval. Then we find some a' > such that Ran(l/(//y P )) C 
D(e a ! x '). In particular, 

|| e a 'l x l l z (//» p ) 11^+^ < const . (5.5) 

J/ G x is modified, then we get uniform lower bounds on a' and uniform upper 
bounds on the constant in (|5.5p . provided that we have uniform upper bounds 
on d_i,di,S np and uniform lower bounds on S np — sup/. 

Proof: The proof is given in the succeeding three subsections. □ 
Note that the potential V is not assumed to be a small form perturbation of 
\J — A in Part (i) of the previous theorem. In particular, the assumptions of 
(i) cover the Coulomb potential — 7/|x| with coupling constants 7 £ [0,7 PF ] 
including the critical one. This improves on [37] where Coulomb potentials 
have been treated, for subcritical 7. By a modification of the arguments of 
this section it is actually also possible to prove exponential localization for the 
no-pair operator with Coulomb potential in the critical case 7 = 7" p , which is 
not covered by Part (ii) of the above theorem; see [2"5] . 

The bound on the decay rate a of Part (i) has been found first in [25] (where 
only the Coulomb potential is treated explicitly). It reduces to the typical 
relativistic decay rate known for eigenvectors of electronic Dirac or square-root 
operators when G x is set equal to zero. 

5.1 A general strategy to prove the localization of spectral 
subspaces 

The general strategy of the proof of Theorem 15.11 is essentially due to [3] . We 
shall present a variant of the argument used in [3] in Lemma l5.2l below. In order 
to apply this lemma to Hy we shall also benefit from some useful observations 
made in |16j . The main advantage of Lemma 15.21 and its earlier variants is 
that it allows to study the localization of spectral subspaces without any a priori 
knowledge on the spectrum. Its proof does not exploit eigenvalue equations as 
it is the case in Agmon type arguments nor does it assume discreteness of the 
spectrum or the presence of spectral gaps. This is important for us since the 
spectra of both the no-pair and the semi-relativistic Pauli-Fierz operators will 
be continuous up to their minima. Roughly speaking the proof of Lemma 15.21 
rests on a combination of the following: 

• The representation (I5.1ip for the spectral projection l/(frf.). Here a com- 
parison operator, Y, enters into the analysis whose resolvent stays bounded 
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after conjugation with suitable exponential weights, for all relevant spec- 
tral parameters. (|5 . 11 1) is valid since its also satisfies li(Y) = 0. 



• The Helffer-Sjostrand formula (re-derived below for the convenience of the 
reader) which is used to represent smoothed versions of Ij(ffi) and li(Y) 
as integrals over resolvents. 

• The second resolvent identity; in fact, Y will be chosen such that — Y 
is well-localized and is hence able to control exponential weights. 

In the somewhat technical parts of this section following after Lemma [5.21 we 
shall verify the applicability of Lemma 15.21 to our models by defining and ana- 
lyzing suitable comparison operators Y . 

Let us now introduce some prerequisites for the proof of Lemma 15.21 In 
order to find a representation of the spectral projection which is accessible for 
an analysis we smooth out the projection and employ the Helffer-Sjostrand 
formula. More precisely, let / C (— oo, S") be some compact interval. Then 
we pick some slightly larger compact interval J C (— oo, S"), J D /, and some 
X G Qf (R, [0, 1]) such that % = 1 on / and x = outside J. We pick another 
cut-off function p £ C^°(R, [0, 1]) such that p = 1 in a neighborhood of and 
p(y) = 0, for \y\ ^ 1/2, and extend x to a compactly supported smooth function 
of the complex plane setting 

X(x + iy) := \(x) + x'( x ) W p{v) , x,y £R; 

compare, e.g., [ITJ[T2]. Then we have 

2dzx(z) := (d x +id y )x(z) = x'(x) (1 - p(y) - y p'(y)) + x"(x) iy p(y) , 

for every z = x + iy € C, and the choice of p implies 

\d*x{z)\ < C x \lmz\, zeC, (5.6) 

for some C x £ (0,oo). Moreover, the following Helffer-Sjostrand formula is 
valid, for every self-adjoint operator, X, on some Hilbert space, 

X (X) = f (X - z)- 1 dp{z) , dfi(z) := --dtxiz) dxdy . (5.7) 

If X were a complex number (|5.7[) were just a special case of Pompciu's formula 
for some path encircling the support of x- X can, however, be inserted in that 
formula by means of the spectral calculus; see, e.g., |12) . 

To facilitate the discussion of operator domains we replace a|x| in (|5.4[) by 
some F : — > R satisfying 

FernlX R ), F = ar °n B r/2 (0), F^O, \VF\^a, (5.8) 

where r 1 is the parameter appearing in (I5.3|) . 
Finally, we recall our notation jj £ {np,PF}. 
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Lemma 5.2 Let I , J, \, and C x be as described above and assume that Y is a 
self-adjoint operator in Jt? with T>{Y) = T>(Hy) and Y > sup J. Furthermore, 
let a > and assume there exist C, C € (0, oo) such that, for all F satisfying 

|| e F {H v - Y) || < C, sup || e F (Y - z)' 1 e~ F || < C . (5.9) 

ze.J+iR 

Then Ran(l/(#^)) C X>(e a|x| ) and 

|| e°W liCHjJO || < CC'C x £( S npp( X '))/7T, (5.10) 

where C denotes the Lebesgue measure on R. 

Proof: Since 1/ = \ 1/ and J C g(Y) we have 

t I (Hl) = (x(H i v )-x(Y))l I (Hl). (5.11) 

Applying the Helffer-Sjostrand formula (|5.7[) and the second resolvent identity 
we infer that 

e F lj(Hl) = f {e F (Y - z)' 1 e- p }{e F {Y - fl* - z)" 1 d^z) . (5.12) 

Jc 

Estimating the norm of these expressions using (|5.6[) and \\{Hy — z) _1 || ^ 
l/|Imz| we obtain (|5TT0|) with a|x| replaced by F. Then (|5TT0|) follows by 
choosing a suitable sequence of functions F n satisfying (|5.8|) and converging 
monotonically on {|x| ^ 2r} to a|x| and applying the monotone convergence 
theorem to e Fn li(Hy) ip, for every ip g Jff. □ 

5.2 Choice of the comparison operator Y 

In the next step we thus have to find a suitable operator Y fulfilling the condi- 
tions of Lemma |5T2"1 The hardest problem is to verify the second bound in (|5.9p 
and this is actually the main new mathematical challenge in the study of the 
exponential localization in our non-local models. We defer the discussion of the 
second bound in (|5.9[) to the next subsection. 

In order to construct an operator Y whose spectrum differs only slightly 
from the spectrum of the free operator Hq, so that Y > sup J, and which is 
defined on the same domain as H v we simply add some bounded term to H v 
which compensates for the singularities and wells of the electrostatic potential 
and thus pushes the spectrum up to the ionization threshold. A similar choice 
of a comparison operator has been employed in the non-relativistic setting in 
[Tfi] . More precisely, we first introduce a scaled partition of unity. That is, we 
pick xo,R,Xi,R e C°°(R 3 , [0,1]), R > r, such that X o,R = 1 on B R (0), Xq ,r = 
on R 3 \ B 2 r(0), xI,r + Xx,r = 1, and ||Vxfc,ij||oo ^ c/R, k = 0,1, for some 
i?-independent constant c £ (0, oo). Then we define Y as follows: 
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The semi-relativistic Pauli-Fierz operator. We define 

Y™:=H™+(X™-E™) X l R , (5.13) 

for some R ^ 1 which shall be fixed sufficiently large later on. Of course, Hy F 
and Yy F have the same domain and the first bound in (|5.9[) which provides the 
control on the exponential weights holds trivially, 

PF y^PFni ^ /v PF R lPF \ll c - F ^2 II s /-v PF j^PF \ ar+2aR 



\e*{H™ -Y™)\\ < (E — Py ) ||e Xo »||oo < (E - £7 



e 



for every P satisfying (|5.8[) . We shall sketch the proof of the condition Y" > sup J 
which follows from the next lemma. 

Lemma 5.3 Y FF ^ E PF — o(P°), R — > oo, in the sense of quadratic forms on 
& , where the little-o-symbol depends only on V and Xo,i- 

Proof: We employ the localization formula Q3.27P with K = 1, the error 
estimate ggg) , and H V F > P PF , H pf ^ E PF , to get 

Yy F ^ Xo,R H FF xq,r + Xi.fl Hq Xi,R 

+ (E PF - P PF ) X ^ + V - 0(1/ R 2 ) > E PF - o(P°) . 

We also used that sup| x |^ H |V(x)| = o(R°) and Xo,/j + Xi,i? = !• 1=1 
P/ie no-pair operator. Since |-Da| = -Pa + -Pa I-^a| we have 

H v p = \D A \+P+VP++ ff f diag , ff f diag := P+ Hi P+ + P A JT £ P A , 

on 0. We write Hy P R := PP^ P - (1/-R) # f diag , for some P > 1, so that P^ := 
inf cr[P^ p fi P+] > -oo by (13311) . and define 

F7 = + (E"p - E» P R ) X o,r Pi Xo,r ■ (5.14) 

Again it is clear that H v p and Yy P are self-adjoint on the same domain. Also 
the first bound in (|5.9p again follows trivially, 

\\e F {H v p - F7)|| < (E n P - E V P R ) e ar+2aIi , 

for every F satisfying (|5.8[) . Besides the second bound in (|5.9p it remains to 
derive the following lemma. 

Lemma 5.4 Y^ p ^ E n P - E n P/P- o(P°) -kdj/R 2 as quadratic forms on 9, 
where k £ (0, oo) and the Little-o-symbol depend only on V and Xo.i- 

PROOF: Again we use an IMS type localization formula to infer that 

Y* p > X o,R H V P R X o,R + Xi.R H»% X i,r + (E np - E^ R ) X o,R P£ Xo,R 

+ £ ^ + Xi,RPl V P+ X i,R + [**.« . lXk,R, Hy P R }] . (5.15) 

fc=0,l 
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As a consequence of (12T251) . ([236]) . (1X2511 . (15111) . (15321) . and # f 2£T f diag the 
double commutator in the last line is bounded from below by -{k/R 2 )(H? lag + 
d\ + 1), where k € (0, oo) depends only on V and Xo.i- To control this error 
we use the term ^Hf lag in (15.151) . Furthermore, we put hr := Xo,i?/2j so that 
Xi,M*R = 0. Then (1 - p£) V = o(Ji°) and, by (£29, 

-Xi,**a ^ Xi.fl < -Xi,fl[^A » ^ [MR, il] + X?,fl 
< + 0(1/ R 2 ) (ff f diag + d 2 1 + l), 

so that the second term in the last line of f|5 . 15[) can again be controlled by 
the first one. Using these remarks, H V P R > E V P R P+ + (1 - 1/i?) £ n P P^, and 
Hqji >0-~ 1 / R ) S " P ( b y we arrive at the assertion. □ 

5.3 Conjugation of Y with exponential weights 

In order to prove Theorem 15.11 it only remains to verify the second bound in 
(15. 9|) . The following lemma (37] gives a criterion for this condition to hold true. 

Lemma 5.5 Let Y be a non-negative operator in J^f which admits & as a form 
core. Set b := inf u{Y) and let J C (—00,6) be some compact interval. Let a G 
(0, 1) and assume that, for all F satisfying ([5.8p . we have e ±F Q(Y) C Q(Y). 
(Notice that e ±F maps into itself.) Assume further that there exist constants 
c(a), f{a), g(a), h(a) G [0, 00) such that c(a) < 1/2 and S := b — max J — b g(a) ~ 
h(a) > and, for all F satisfying (|5.8[) and ip G , 

\{<p\(e F Ye- F -Y)<p)\ < c(a) (tp\Y ip) + f(a)\\<pf, (5.16) 

Re{tp\e F Ye- F tp) ^ (1 - g{a)) ( ip \ Y ip ) - h(a) \\tp\\ 2 . (5.17) 
Then we have, for all F satisfying (|5.8|) , 

sup \\e F (Y-z)- 1 e- F \\ < S' 1 . (5.18) 

ze.J+iR 

Proof: We only sketch the proof and refer to Lemma 5.2 of [37] for the 
details. The assumptions e ±F Q{Y) C Q(Y) and (j5l6l) ensure that the closure, 
Yp, of (e F Ye~ F ) \@ agrees with the closed operator e F Ye~ F . The bound 
f|5 . 1 T[) shows that the numerical range of Yp is contained in the half space 
{z G C : Re 2: ^ sup J + 5}. Moreover, we can argue that, for z G J + iR, the 
deficiency of Yp — z is zero and, hence, the norm of (Yp — z)^ 1 = e F (Y — z)e~ F 
can be estimated by one over the distance of z to the numerical range of Yp. □ 

The semi-relativistic Pauli-Fierz operator. Next, we apply Lemma 15.51 to 
Hy F . In order to verify condition (|5.17j) with a good bound on the exponential 
decay rates we apply the following technical lemma from [25] : 

Lemma 5.6 For all a G (0, 1), F satisfying (j5.8|) . and ip G @>, 

Re{ip\e F \D A \ e~ F V > > < <p \ (D\ - | VF| 2 ) 1/2 ip ) 

>{<p\(\D A \-p(a))<p). (5.19) 
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Proof: For every ip G @, we infer from (|4.12|) that 

< ? | ( e F |23 A | e" F - (^i - |vf| 2 ) 1/2 ) v ) = / jy- r,] ?-*L , 

Jo 7r 

with 

Jfo>; rj\ :=(<p\ (M F {r,) - e F ^ («) e~ F ) <p ) , 
&g(v)--=(D 2 a -\VG\ 2 + V )-\ Ge{0,F}. 

Now, let 4> '■= e F (D A + rf) e~ F ip, for some ip G ^. Then 

Re ( | e F ^ (r?) e" F <?(>) = Re < e F (D 2 A + 77) e~ F V> | </> ) 
= ((^I-|V^| 2 +r/)V|^) 
^ {l-a 2 +rj)U\\ 2 

Since is essentially self-adjoint on 3! and multiplication with e~ F maps ^ 
bijectively onto itself, we know that (D A + rf) e~ F S> is dense in 34?. Since F is 
bounded we conclude that the previous estimates hold, for all <fi in some dense 
domain, whence Re \e F e~ F ] as a quadratic form on Next, we set 
Q := (a ■ VF) D A + D A (a- VF) and let 

<P ■= (D 2 A - |VF| 2 + rf)iP = e ±F (D 2 A + rf) e^ F ^ T iQ^P, 

for ip G ^. Then 

Jfo n] = *( e'^oiv) e F y> | Q V ) = *< V- 1 Q 1> ) + < Q V> \ e F ^o(v) e~ F Q ip ) . 

Here D A — \WF\ 2 is essentially self-adjoint on ^ and Q is symmetric on the same 
domain. Hence, Re J\tp\rf\ ^ 0, for all <p in a dense set, thus for all (p G J^, 
and we arrive at the first inequality in (|5.19|) . Since the square root is operator 
monotone, |VF| $S a, and \D A \ ^ 1, we further have 

(Dl - |VF I 2 ) 1 / 2 ^ |£> A | + (Ui - a 2 ) 1 ! 2 - \D A \ > |D A | - p(a) • 

□ 

In what follows we abbreviate 

K F := [F+ , e F ] e" F , 

and recall from (|2~24| that || |L» A | 1/2 £f || = O ao (a), for all F satisfying (JSHJ 
with < a ^ ao < 1. 

Lemma 5.7 For aZZ < a ^ ao < 1 and F satisfying (|5.8[) . 

Rc(^|e F yP F e - F ^)^(^|rP F ^)-p(a)||^|| 2 , ^G^. (5.20) 

Moreover, for all e > and ao G (0,1), i/iere is some constant, C(ao,£, V) G 
(0,oo), sztc/i i/iai, /or a/Z F satisfying (|5.8|) wii/i a G [0, ao] and <p G @, 

|(^|(e F yP F e - F -Y* F )<p)\ ^ e(v\Y$ F ^) + C(a ,e,V)M 2 . (5.21) 
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Proof: (|5.20j) follows immediately from (|5.19|) . To derive (|5.2ip we write 

e F \D A \ e- F -\D A \ = -2 D A JC F + ia ■ (VF) e F S A e" F 

on 3 and make a little observation. Since F = ar on B r / 2 (0) we find some 
jti € C°°(R^ [o, 1]) such that /i = 0on B r/4 (0) and VF = fi VF. For 99,^ € 3>, 
we thus have (recall (j2~18f and ([2~22|» 

|(£) A¥3 |/C F V>| «S / \{^\D A R A {iy)^ia-VFR F A {iy)i,)\^-. 

Here we can write 

D A R A {iy)^ = fi\D A \ S A R A {iy)+[D A , fi] R A {iy)+D A R A {iy) \p,D A ] R A (iy) , 

where ||[-Da,/^]|| ^5 ||V/i||oo = O(l)- On account of 

\\\D A \V 2 R A (iy)\\ ^(l)(l+y 2 r 1/4 

and 1 1 -Da -Ra(*2/)|| = 0(1) it is now straightforward to verify that 

\(D A ip\K F ^)\ < O ao (a) {|| |Fa| 1/ V^|| + IMI} W- 

Some elementary estimates using HVFHoo ^ a, ||e F SAe~ F | = O ao (l), and the 
previous bound now show that 

|(^|( e F |£>A|e- F -|£>A|)^>| 

< ei < ^ V | \D A \ + {e^ 1 O ao (a 2 ) + O ao {a)) \\v\\ 2 (5.22) 

< si 0(1) (n<p\ Y FF flip) + const(a ,£i) \\p\\ 2 , 

for every e\ € (0,1]. In the second step we used that \x V is bounded because 
fi = on B r /4(0). Since we may assume that there is some jus S C^°(IR 3 , [0, 1]) 
such that /x 2 + Jl 2 = 1 we can employ an IMS localization formula as in the 
proof of Lemma [531 to show that /i Y FF [i Fy F + 0(1) on f^. Altogether this 
proves (|5~2Tj) . □ 

Lemma 5.8 There exist constants, c\-,c% € (0, oo), such that, for all a S 
(0, 1/2] and ±F satisfying (I5.8[) . and p E 3 1 , 

(e F p\Y FF e F p) < c 1 || e F || 2 (^|rP F ^)+ C2 ||e F || 2 ||^|| 2 . (5.23) 

In particular, e F Q(Y FF ) C Q(Y^ F ). 

Proof: We pick a smooth partition of unity with respect to the electron 
coordinates, /i 2 , + /i 2 = 1, where supp(/io) C i?r/2 (0) and /j,q = 1 on i? r /4(0). 
Then ( e F p \ Y^ F e F <p ) = ^=o,i ( Mi e ^ V | ^i/ F Mi e F <p) + R v , where | < 
0(l)||e F || 2 II^H 2 . (This holds in particular for F = 0, of course.) Therefore, 
it is sufficient to prove the bound (I5.23|) with p = m ip, i — 0, 1, ip £ < 2>. For 
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p = fiQ ip, the bound holds, however, true trivially, for all c±, c 2 > 1, since F = 1 
on the support of /lo- 

Let us assume that p = fii ip, for some ip € @, in the rest of this proof. 
Of course, \\xo,R e F 'fill 2 ^ ll eF || 2 \\Xo,R ^ll 2 and, since Hi and e F commute, 
\\Hl /2 e F p\\ 2 < ||e F || 2 ||.ff f 1/2 ^|| 2 . Furthermore, K^|V^)I < 0{l)\\p\\ 2 , since 
V is bounded on supp(/ii). To conclude we write \Dx\ = P£ Da — P A -^A and 
employ the following bound derived in |39[ Equation (4.24) and the succeeding 
paragraphs] , 

(p\e F P±(±D A )e F p) ^c 3 \\e F \\ 2 (p\P±(±D A )p)+c 4 \\e F \\ 2 \\p\\ 2 , (5.24) 

for every tp £ < 2>. We actually derived this bound in 39 for classical vector 
potentials. The proof works, however, also for the quantized vector potential 
without any change. Moreover, we only treated the choice of the plus sign in 
(|5.24[) . But again an obvious modification of the proof in [39] shows that (|5.24[) 
is still valid when we choose the minus sign. □ 

The no-pair operator. The following lemma implies that the conditions (|5.16p 
and (|5.17[) are fulfilled in the case of the no-pair operator, too. 

Lemma 5.9 There exist c = c(V) € (0, oo) and c' = d(V, d_i, d 1: T, np -E V P R ) G 
(0,oo) such that, for all F satisfying (|5.8p . 

\(<p\(e F Y^e- F -Y»*)<p)\ < 0(a) (<p\ {cY? + J) p) , <p£®. 

Proof: On account of (|5.22j) and \\e F X o,r P a Xo,r e~ F - Xo.R P A Xo.rW < 
\\ICf\\ = O ao (a) it suffices to consider 

A ± (T) := e F P^T P± e~ F - P^ T P± = 2Rc [P± T SP] +SPT SP, 

where SP := e F P^ e" F - P^ and T is H f or V. Clearly, 

\(p\A ± (T)p)\ < e(p\P^TP^p) + (l + e- 1 )\\\T\ 1 / 2 SP^\\ 2 , 

for all e > and ip e 9. Since, by (|235l) and (j236l) . || \T\ X I 2 SP pf sC 
0{a 2 )(p\(H { + (4di) 2 + l)(p) ^ 0{a 2 )(ip\ (ff f diag + d\ + 1) p ), we may choose 
e proportional to a and use p. 311) to conclude. □ 

To complete the proof of Theorem 15 . 1 1 also in the case of the no-pair operator we 
note that the bound (|5.23p still holds true when Yy F is replaced by yg np . To this 
end we only have to observe in addition to the remarks in the proof of Lemma lSTSl 
that \\Hf /2 P^e F p\\ ^ 0(1) ||e F || ||i? f 1/2 p\\. This follows, however, immedi- 
ately from (HU) which implies ||iJ f 1/2 P± e F p\\ ^ (1 + ||C 1/2 ||/2) \\e F tl\ 12 p\\. 
Thus, e ±F Q(Yq P ) C Q(Yq P ). By Theorem 13.61 and the assumptions on V we 
know that Q(Yy P ) = Q(Yq P ) and we conclude. 
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6 Existence of ground states with mass 



In this section we present an intermediate step of the proof of the existence of 
ground states for H v . Namely, we prove the existence of ground state eigen- 
vectors, (j)^, for modified Hamiltonians, H Vm , which are defined by means of 
an infra-red cut-off coupling function. The infra-red cut-off parameter, m > 0, 
is referred to as the photon mass. Later on in Section [8J we shall remove the 
infra-red cut-off by showing that every sequence, {<^i.}j, rrij \ 0, contains 
a strongly convergent subsequence whose limit turns out to be a ground state 
eigenvector of H v . The compactness argument used to show this in Section [8J 
requires the infra-red bounds derived before in Section [7] 

In the present section the existence of $ m is shown by discretizing the photon 
degrees of freedom. After the infra-red cut-off operators Hy m have been defined 

in Subsection 16.11 we construct discretized versions of them, denoted by H Vm , 
in Subsection 16.21 We collect some technical estimates needed to compare the 
original, infra-red cut-off, and discretized operators in Subsections 16.31 and 16.41 
As another preparation we study the continuity of the ground state energy and 
ionization threshold with respect to the parameters m and e in Subsection 16.51 
The main result of this section, Theorem 16. 1 II on the existence of 4>l, i: is stated 
and proved in Subsection l6 . 61 and we refer the reader to that subsection for some 
brief remarks on its proof. Many arguments of this section (in particular those 
in Subsections I6.4l and l6.5|) are alternatives to the corresponding ones in [2"5I 125] . 

In the whole section G x is the coupling function given by (I2.2|) . To clarify 
which properties of V are exploited we introduce the following hypothesis. It is 
fulfilled by the Coulomb potential in the subcritical cases: 

Hypothesis 6.1 In the case jj = PF the potential V <E Lf oc (R 3 , R) is relatively 
form-bounded with respect to \/—A with relative form bound strictly less than 
one. In the case (j = np the potential V € Lf oc (R 3 , R) satisfies H 1 (R 3 ) C T>(V) 
and (|3.30j) with a < 1. 

We shall strengthen the assumptions on V later on in order to apply the 
localization estimates of Section [5] 

6.1 Operators with photon mass 

For every m > 0, the infra-red cut-off coupling function is given as 

G x , m (fc) := -e 1{? "^ A} e -*' x s(k) , (6.1) 
27ry |k| 

for all x £ R 3 and almost every k = (k, A) G R 3 x Z2. To compare G XyTn with 
G x defined in fl2.2j) we introduce the parameter 

A 2 (m) := f (u(k) +u(k)- 1 ) sup|G x (fc) - G^ m {k)\ 2 dk . (6.2) 
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Of course, A 2 (m) = (e 2 /27r 2 ) Jj k , <m (l + |k|- 2 )d 3 k -> 0, as m \ 0. For m > 0, 
we further define the infra-red cut-off vector potential, 

A m :=a t (G m ) + o(G m ), a B (G m ) := / 1 C 4 <g> a'(G x , ra ) d 3 x , 

and the infra-red regularized Hamiltonians 

H^ F m := \D A J + V + H { , (6.3) 

flSfm : = P X„ ( D A m +V + Hf) P+ m . (6.4) 

We define these operators as self-adjoint Friedrichs extensions starting from Ql. 
The ground state energies and ionization thresholds, for positive photon mass 
m > 0, are denoted by 

E\ m := inf a[H^J , £»„ := inf <r[H* >m ] , 

As a first step we introduce a truncated Fock space where the radiation field 
energy Hf is bounded from below by m > on the orthogonal complement 
of the vacuum sector. Namely, we split the one-photon Hilbert space into two 
mutually orthogonal subspaces 

jr = Jtr>®jr<, J(T> := L 2 (A m x Z 2 ) , An := {|k| > m} . 

It is well-known that J^JT] = ^b[Jf>] ® &b[Jf<]. We observe that A m 
creates and annihilates photon states in J(f> only and iff leaves the Fock space 
factors associated with the subspaces Jifm invariant. We shall designate oper- 
ators acting in the Fock space factors or by a superscript > 
or <, respectively. Under the isomorphism 

& = (L 2 {^x A )®^ h [j(r>])®& h [je<] =■. m>> ®& h [j*r<], (6.5) 

we then have D Am = D A > ® 1, \D Am \ = |-D A >I ® 1 , p Am - P a> ® 1 ' and 
i? f = flf + ff< with Hj m := dr(u\ AmXZ2 ), Hf m := dT (wf^zj. As a 
consequence, the semi-relativistic Pauli-Fierz and no-pair operators decompose 
under the isomorphism (|6.5I) as 

= ff^o i 1 + 1 i H< , (6.6) 

= H% mfi ®l + P+>®H<, 

where 

#S,m,o := \D A >J+V + H>, 

Hy P m , := P+> (£> A > + V + H>)P+> . 

The latter operators act in the Hilbert spaces Jt?^ and P^> , respectively. 
The following lemma [28, 29] shows in particular that it suffices to prove the 
existence of ground states in these truncated Hilbert spaces. 
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Lemma 6.2 Assume that V fulfills Hypothesis \ 6. 11 Then, for all m > 0, 
4,™ = inM< m , ] , 4, = inf <r[< m , ] • 

Moreover, if E Vm is an eigenvalue of H Vm0 , then it is an eigenvalue of Hy , 
too. 

Proof: It is clear that H Vm0 ® 1 < H Vm , whence mfa[H v ] < 

infer [if y ]. Next, we pick a minimizing sequence of normalized vectors ip£ £ 

Q(Hv, m ,o), (^n\ H V,m,o^n) -> infa[^ m0 ]. Setting tf, n := ^> <g> fi>, where 
Q > is the vacuum vector in j?k[j£^], we observe that 

(lPn\Hy^ n ) = (lP>\Hl m ^>), 

thus infcr[ffy m0 ] ^ inf ir[-ff|> m ]. Likewise, if $ m £ M£ is a ground state 
eigenvector of H Vm , then 0^ ® fP is a ground state eigenvector of H Vm . □ 

In order to show the existence of a ground state for H v in the next step it is 

sufficient to show that the spectrum of H Vm Q is discrete in a neighborhood of 

E Vm . [Ey m is contained in the essential spectrum of H Vm on the contrary.) 
A general strategy to achieve this would be the following. We could seek for 
a self-adjoint operator, A, satisfying — oo < A ^ H Vm0 and having discrete 

spectrum in (— oo, Ey + c], for some c > 0. If such an operator A exists 

then also H Vm has discrete spectrum in (— oo, E Vm + c]. We need, however, 

a modification of this strategy. Let x denote the spectral projection of H Vm Q 

corresponding to some half-line (-co, E Vm + c], c > 0. Then we seek for a 
self-adjoint auxiliary operator A such that 

xAx^X {< m , ~ E V m -2c}x^~c X (6.7) 

and Trjx^lx} > — oo, where Tr denotes the trace. For in this case we have 
Tr{x} < oo. The latter strategy is advantageous since we only have to compare 
A and H v m on the range of x whose elements are exponentially localized by 
Theorem 15.11 provided that c > is appropriately chosen. A suitable com- 
parison operator A is constructed by means of a discretization of the photon 
momenta in the next subsection. 

6.2 Discretization of the photon momenta 

On we introduce a discretization in the photon momenta: For every e > 0, 
we decompose A m = {|k| ^ m} as 

An= (J Q S m {v) , Q'mH ■= (V +{-6/2, E/2) 3 )nA m , V£ ( £ Z) 3 . 
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Of course, for every k G A m , we find a unique vector, £ E (k) G (^Z) 3 , such that 
k G Q^(i/ E (k)). To each v G (eZ) 3 with Q £ m (v) /0we further associate some 

Xm,e( U ) ^ Qmiv) SUCn tnat 

km, £ (^)| = inf |k| . 

In this way we obtain a map 

t'e : An x Z 2 — > R 3 , fc = (k,A) .— ► i/ e (fc) := x m , e (P e (k)) . (6.8) 

It is evident that the vectors x r n,e{ l/ ) can be chosen such that 

v e {— k, A) = — i/ e (k, A) , for almost every k G A rn ■ (6-9) 

The set of Lebesgue measure zero where the identity (|6.9p might not hold is 
contained in the union of all planes which are perpendicular to some coordinate 
axis and contain points of the lattice (e Z) 3 . We define the e-average of a locally 
integrable function, /, on A m x Z2 by 

[Psfm := I /(p ' A)d3p ' (fU0) 

and introduce the following discretized coupling function, 

G„(k) := -(e/2n)e-^ k > x P^Ia,^- 1 / 2 e](k) , 

for all x G R 3 and almost every k = (k, A) G A m x Z 2 . In order to compare 
G x ,m with G Xim)£ we put 

A 2 (a,m,e) := / + oT 1 ^)) sup {e^W |G x . m (fc) - G x ,„^(fc)| 2 }dfc, 

-A m xZ 2 

(6.11) 

for a, m, e > 0. It is elementary to verify that A* (a, m,e) — > 0, as £ \ 0, for 
fixed a, to > 0. Notice that A*(a,TO,e) did not converge to zero if we chose 
a = 0. The fact that we need some weight function in x to control the difference 
between G x . m and G Xim . e is one of the reasons why a localization estimate 
for spectral subspaces is required to prove the existence of ground states. The 
discretized vector potential is now given as 

A m , e :=a^(G mte ) + a(G m , E ), a J (G m , £ ) := / 1 C 4 ® a»(G w ) d 3 x . 

Jr 3 

The reason why we choose vectors v e fulfilling (16.91) is its consequence 

Gx,m : e(~ k, A) = G X;rnj£ (k, A) . 

The latter identity ensures that different components, Am| E (x), A^fy), x, y G 
R 3 j hj G {1,2,3}, of the discretized vector potential still commute. We have 
used this property in Section [3] The dispersion relation is discretized as 

tj e (k) := inf { |p| : p G Q^(£ £ (k)) } , k — (k, A) G A m X Z 2 . 
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Then \u s \ ^ lu £ on A m x Z2 and 

to <w e w on ^ m x Z 2 , flf, m ,e : = dT(u! £ ) < ff f ^ m . (6-12) 

Here the operators Pf, m ,e and Pf^ are acting in J?b[J^]. Finally, we define 
discretized Hamiltonians, Py m e , acting in L 2 (R 2 , C 4 ) ® ^[JO^], 

H v F m ,e : = I^A m , E I + V + Pf, m>£ , 

^y P m , £ := ^A m , e (^A m . E + V + P f>m , E ) P+_ . 

6.3 Comparison of operators with different coupling func- 
tions 

In order to compare the various modified operators we derive some bounds on 
differences of projections whose proofs are essentially consequences of the ideas 
collected in Subsection 12 .4! and the bounds 



|| a ■ (A - A m ) H f 1/2 || = C(A(m)) , (6.13) 
||a-(A> -A m , e )i? f ~ 1/2 e- a l x l || = 0(A*(a,m,e)). (6.14) 

Here we use the notation (|6.2[) and (|6.1ip . 

Lemma 6.3 Let V be a symmetric multiplication operator in i 2 (R 3 ) which is 
relatively form bounded with respect to y— A. 

(%) Sei i/f = iff + E, for some sufficiently large E ^ 1 depending on e and A, 
and let v ^ 0. Then, as m \, 0, 

||p A r /2 (Pi-Pi m )^ f - 1/2 || < 0(A(m)), (6.15) 

|| I^a„J 1/2 (Pi - Pt„)Pr 1/2 II < 0(A(m)) , (6.16) 

||Pr(PA-PA,J^ 1/2 H < <?(AM), (6.17) 

IH^CPi-P^)^- 1 !! < 0(A(m)). (6.18) 

(ii) Let Hf be H> m +E or H { ^ e +E, for some sufficiently large E 1 depending 
on e and A, and let v and ao G (0, 1) . Then, for every a G (0, ao] and 
P G C°°(IR|, [0,oo)) satisfying |VP| sC a, P(x) ^ a|x|, /or a// x G R 3 ; and 
P(x) = a|x|, /or /arc/e |x|, and for all sufficiently small m,e > 0, 

I|Pa>| 1/2 (^a> - p l m , e )^f" 1/2e " F H < 0(A,( a ,m,e)), (6.19) 

||I^A m ,J 1/2 (P^> -P± m J^ f - 1/2 e- F || ^ 0(A,(a,m,e)), (6.20) 

ll^f'( P A> - p A m<e )^ 1/2 e- F || < 0(A,(o,m,e)), (6.21) 

|| \V\ 1/2 (Pj£> - Pj m J Pf" 1 e" F || < 0(A,( a ,m,e)). (6.22) 
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PROOF: By the assumption on V and Theorem l3.4l we have |V| < C (|Da| + 
H f ). Therefore, the bounds (f6TT8l) and (|Q2|) are consequences of ([6715]) fc ((F!T7|i 
and (|6.19[) fe (|6.21l) . respectively. In order to prove the remaining estimates we 
pick two vector potentials, Ai and A 2 , such that the set {Ai, A 2 } equals either 
{A, A m } (in which case F := in what follows) or {A> , A TOi£ }. For j = 1,2 and 
y e R, we set R A] (iy) := (D Aj - iy)- 1 and (iy) := (D Aj +ia-VF- iy)- 1 . 
(Recall (|2.20l) and (|2.2ip .) Then we have the following resolvent identity, for 
(j, ^ 1/2 and ye R, 

(R Al (iy)-R A2 (iy))e- F Hf» 

= R Al (iy) a • (A 2 - A x ) Hfe~ F Rl 2 (iy) T F (iy) , 

where T F (iy) is the bounded operator defined in (|9.5[) below (with A replaced 
by Ai). Using (I2.18|) we find, for all 4> and / in the Hilbert space, 

\{f\{Pt-Pi 2 )H^e- F <i>)\ (6.23) 

< f \(f\R Al (w)a-(A2-A 1 )H-^e- F Ri 2 ( l y)T F (zy) < p)\^-. 

In the case fj, = 1/2, we choose / = \D Al \ 1/2 ip, ip € V(\D Al \^ 2 ), and observe 
that the integrand in (I6.23[) is then bounded by 

0* WIDa^RaMW \\R A2 , L (iy)\\ ||T 1/2jL (%)|| U\\ uw, 

which is integrable due to (|2.21j) . (|6. 13[) . (|6.14l) . and the spectral theorem. Here 
and below 0* = 0(A(m)) or 0* = 0(A* (a, m, e)) depending on the choice of 
Aj. This concludes the proof of (|6T5|l . (|6l6| . f[67T9]) . and (|6~20| . 

In order to prove (|6.17|) and (|6 . 2 1 1) we infer from Lemma l9~Tl that the com- 
mutator T v :— Hf [a ■ (A 2 — Ai) e~ F , Hf v ] extends to a bounded operator 
with ||T„|| = 0*. Together with ([6TT5]) , ([STTi|) . and (JUTJ) this implies, for i/ > 
and ip,(f> £ V(H?), 

| ( | i? Al (iy) o • (A a - Ai) H; u ~ 1/2 e- F R F A , (iy) T F +1/2 (iy) ) | 
= \{i>\R Al (iy)r F (iy)H?OL-(A. 2 -A l )H- v - l, ' 2 e- F x 

xRl 2 (iy)r F +1/2 (iy)<j>)\ 
^ 0, Pa, Mil ||Tf ||< 2 ^)|| ||T^ 1/a (iy)|| ||V|| U\\ 
^O^l + y 2 )- 1 U\\H\\. 

Therefore, (|6 . 1T[) and (|6 . 2 1 1) follow from ()6.23|) upon choosing fi = v + 1/2 and 

f = H?il>. □ 

6.4 Higher order estimates and their consequences 

As a preparation for the proof of the existence of ground states for H v we 
derive bounds on certain expectation values of the difference H v m0 — H v in 
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this subsection. Moreover, we compare the ground state energies and ionization 
thresholds of Hy m with those of H v m e defined by 

for 771, e > 0. It is not possible to compare the no-pair operators H v p m and 
Hy m e in a quadratic form sense. For some error terms in the difference of 

these two operators "have the size of |x| i?^ 2 " and can only be controlled when 
we take expectations with respect to states in some low-lying exponentially 

3/2 

localized spectral subspace. To control the higher power H f ' of the radiation 
field energy we need, however, yet another non-trivial ingredient, namely the 
higher order estimates of the next theorem. Because of lack of space we cannot 
comment on the proof of Theorem l6.4l and refer to [5S] instead. We remark that, 
for the semi-relativistic Pauli-Fierz operator, higher order estimates have been 
obtained earlier in |14j . Their proof given in |36) is, however, different and more 
model-independent so that the no-pair operator can also be treated along the 
same lines in [35] . In the case of the no-pair operator only the Coulomb potential 
is considered in [35]. An inspection of the proofs in [35J shows, however, that 
they immediately extend to all potentials satisfying Hypothesis 16.11 

Theorem 6.4 Let e £ R, A £ (0, oo), and assume that V fulfills Hvvothesis \6.1[ 
Then Pffff^J"' 2 ) C T>(H^ 2 ), for every n £ N, and there exist constants, 
£o, mo, C £ (0, oo) 7 such that, for all e £ [0, £o] and m £ (0, mo], 

KZj Plm ^> (H v P miS ~ (E v % t£ - l)PL,V n/2 \\ < W + l^,J) 2 "> 

\\ H f,m,e( H vfm,E - i E V\m,e ~ l )) '^H ^ C ( l + l^m^l) 2 "- 

IfV — 0, then E Vme has to be replaced by £ in these bounds. Analogous 
bounds hold for Hy . 

Lemma 6.5 Let e £ R, A £ (0, oo), and assume that V fulfills Hypothesis 1 6. 11 
Then we find some mo > such that the following holds: 

(i) For all m £ (0, mo] and ip + £ Ran(l(_ co .stt+i] (H v )), 

\(^ + \< m ^ + )-(^ + \HU + )\ < const(SM4|) (m°)||^+|| 2 . 

(ii) For all m £ (0, mo] and il> + £ Ran(l,_ oo s » i+1 i (H Vm )), 

\(i> + \K™^ + )-(^ + \ H U + )\ < const(SL,l4, m l)«(™°)ll^ + l| 2 . 

Proof: We treat only the no-pair operator explicitly. On account of the 
formula | Da | = P£ Da — Pa Da it will then be clear how to obtain the result 
also for the semi-relativistic Pauli-Fierz operator. We remark only once that, 
for instance, the inclusion Pj£ Ran(l(_ oc (Hy)) C Q(H Vm ) follows from 
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the characterization of the form domains in Theorems 13.41 and 13.61 In the rest 
of this section we shall use similar remarks without further notice to simplify 
the exposition. Let 5P := P£ — Pj^ and M := {-Da, V, iff}. Then we have 

H7 ~ Hy P m = PL a -( A ~ A ™) P X„ + E 1 2Re t P A T SP] -SPTSP} 

TEM 

(6.24) 

in the sense of quadratic forms on Q(|D |) H Q(Hf). Now, let 

i,+ e Ran(l ( _ oc , ss+1] (^)). 

Prom (jnHSl), (EH31), ([6TT5]l -([6TI5 1) , and (joTMl) we readily infer that 

\(^+\(H v p -H^J^+)\ < 0(A(m)) {|| |Ca| 1/2 ^ + || 2 + H^fV'+ll 2 }- 

Here || 1-DaI 1 ^ 2 V^H ^ C(1)IIV ,+ II since -0 + belongs to the spectral subspace 
Ran(l(_oo,sJ+i] (Hy)) and the term containing V is a small form perturbation of 
H v by Hypothesis 16. H and Theorem l3.4l or Theorem l3.61 respectively. Moreover, 
\\H{ if} + \\ ^ const(S, \Ev\) because of the higher order estimates. We can 

argue analogously if ip + belongs to a spectral subspace of H Vm . In this case 
the right hand side of the estimate depends on E m and |£V i?n | since we apply 
higher order estimates for H Vm . □ 

We recall once more that some x-dependent weight is required to control the 
difference between and A m e . If we consider only vectors ?/> + in a spectral 
subspace corresponding to sufficiently low energies of H Vm e , then we can bor- 
row this weight from the exponential localization of ip + . At this point we have 
to introduce further conditions on the potential in order to guarantee that there 
are non-trivial spectral subspaces below the ionization threshold. 

Hypothesis 6.6 V satisfies (|5.3j) and there exist c,m*, > such that, for 
all m € (0, m*] and all e € (0, e+] , 

£»-4^c, E«,-4 im >c, Zl tS -E v<mie Zc. (6.25) 

Examples of potentials fulfilling the previous hypothesis have been found in 
Theorem 14.11 In fact, as already pointed out there, the proof of Theorem 14.11 
works also for the discretized operators when the unitary transformation em- 
ployed in Section |4] is modified suitably; see [2H1 [29] for details. There are, 
however, potentials fulfilling Hypotheses 16.11 and 16.61 which are not covered by 
Theorem 14. 1[ for instance, those mentioned in Remark 14.21 This is the reason 
why we work with the very implicit Hypothesis 16.61 in what follows. 

Lemma 6.7 Let e 6 R, A £ (0, oo), and assume that V fulfills Huvotheses \6.1\ 
and \ 6.6\ Then there exist £o,too > such that the following holds, for all 
e G (0, eg] and m £ (0, mo] .' 
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(i) For all A € (£«, m ,£y and ^+ € Ran(l ( _ OC:A] (H^J), 

\{^ + \ H V,m^ + )-{i> + \ H V, m ^ + )\ 

< const(sL, l4,ml. ( E l - ll^ll'- (6-26) 

(ii) For all A € (4 >ro>£ , E?^) and V+ € Ran^.^Ctf^J), 

\(i' + \Hl m ^ + )-(^\Hl nh0 ^ + )\ 

< const(£^ £ , |4„^|, (5t i£ - A)" 1 ) ( e °) ||^+|| 2 . (6.27) 

Proof: Again we only treat the no-pair operator since the proofs for the 
semi-relativistic Pauli-Fierz operator will then be obvious. 

(i): We have a formula for Hy P mfi -Hy P m £ similar to (joTMl) with (A, A m , H { ) 
replaced by (A> , A TOi£ , Hf" m ) where one additional term has to be added, 
namely P+ m e dT(ui - ui E ) Pj^ s - Using this formula, (|oTT4l) . (EHJ)- (EH, and 
| a; — u; £ | ^ y/3e ^ (\/3e/m) uj e ^ (V3e/m) w, for |k| ^ m, which yields 

\(ij + \p+ m jr(u-u e )p+ m _jj+)\ < o^Hi^Vf, 

we arrive at 

\(^\(H^ m ,o-K p m ,Ji' + )\ 

< (0(A,(a,m,e)) + o(e )) {|| |£ A > | 1/2 ^+|| 2 + ||e F H f V+|| 2 }, 

for -0 + g Ran(l(_ 00 X ] (H Vm ))- Here we further have 

2\\e F H i ^\\^\\e 2F ^f + \\Hfyj + \\ 2 , 

where the norms on the right side can be controlled by our exponential localiza- 
tion and higher order estimates, respectively. Part (ii) is derived analogously. 
The dependence on A, the ionization thresholds, and the ground state energies 
of the constants on the right hand sides of (|6.26p and (|6.27p stems from the 
constants in the exponential localization and higher order estimates. □ 

6.5 Continuity of the ionization thresholds and ground 
state energies 

To make use of the bounds of Lemmata 16. 51 and 16. 71 we still have to verify that 
the functions m EjJ^, m i-» E Vm , and e ^ S^ E , s i-» Ey mE are (semi- 
continuous at 0. The continuity of e i-> E VmE will also enter more directly 
into the proof of the existence of ground states in the next subsection. 

Corollary 6.8 Assume that V fulfills Hypothesis \6.1\ Then it follows that 
liirim-^o Ev x m = E v and lim ™^o Sf n = S". 
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Proof: Using Lemma H3Ji) it is not difficult to derive the bounds Ej^ ^ 
E» + o(m°). (Given e > we pick some ip+ £ Ran(l [stl . s « +£) (H^)), = 1, 

and plug it into the quadratic form of Hq m . In the case of the no- pair operator 
we also have to observe that \\P A — > 1, to \ 0.) Since this gives an upper 
bound on E^„ which is uniform, for small m, we can then control the constants 
on the right hand side of the estimate in Lemma [6T5T ii) to get E* ^ E^ + o(m°) 
by a similar argument. Since the results of Section [3J provide uniform lower 
bounds on Ev, m we can now employ Lemma 16.51 in a similar fashion to show 
that lim„ w0 E^ m = E\. □ 

In order to compare the ionization thresholds Y$ m and Y$ m e we need a differ- 
ent argument since x-dependent weights are required to control the difference 
between A> and A m e but the spectral subspaces of the free operators H\ m e 

are not localized. Here the essential self-adjointness of Hq e asserted in The- 
orem [317] is helpful. To be able to work in one fixed Hilbert space we set 

Km,e ■= HZn, E + ^A m , e (I^A m . e | + H {<m , e ) . (6.28) 

Then Lemma IQ1 below implies that E£P £ = inf <j[Hq P 7JI e ] = inf a[HQ P m e ]. 

Lemma 6.9 e — > Hf^ and s — > ffp^ in i/ie strong resolvent 

sense, as e \ 0. In particular, limsup^o s ^ ^L- 

Proof: Since all involved operators are essentially self-adjoint on 9l it 
suffices to show that H%J ne p H^ j0 (p and H^ p m e p ->■ HQ P m0 p, for every 
fixed <p € S>. Since ip has only finitely many non-vanishing Fock space compo- 
nents and the latter are compactly supported it is clear that flf, m , s p — > H{* pi. 
Furthermore, we write 

\Da> I - \DA m , E | = S A > a ■ (A> - A m , E ) + (5 A > - S Am jD Am ^ , 

where S :— D |-D| _1 denotes the sign function, which permits to get 

\\\D K >\p~\D Am Jp\\ < \\e- F a -(A>-A m ^H- 1/2 \\\\e F Hl /2 p\\ 
+ \\(S A > - S Am jH- 1/2 e- F \\ {\\e F H} /2 D p\\ + ||P f 1/2 a • A m , e e F ^||}, 

where Hf := Hf" + E and E = E(e,A) is sufficiently large, independently of 
e. Using (|6.14p . (|6.19[) . and Lemma l9Tl it is now easy to see that Hq F ^ £ ip — > 

PF 

0,m,0 

To show that also H Q lp m E 93 — > ffj p m ^ it remains to observe that 

II p a> ff f p a> <p ~ p L,e H t P L,. 4 

< II ^> - P a_ ) ^ F || \\e F H? /2 P± >m Hf 3/2 e- F \\ \\e F H? /2 p\\ 

+ \\H t (P±> -P± m jH- S/2 e- F \\ \\e F H? /2 p>\\ 0, (6.29) 



P, PF 



44 



and 



(6.30) 

In fact, (j6~29t is valid because of tfttty . (JOJ), and since e F H^ 2 P^ > Hf 3/2 e~ F 
is well-defined and bounded as a consequence of (I2.24[) and Lemma 19. 31 More- 
over, (|6.30[) holds true since \u> — u> E \ ^ \f~Ze and the photon momenta are m 
in modulus on Jt^, and since \\Hf~ -Hf _1 || is bounded uniformly in e > 

by Lemma [231 n 

Corollary 6.10 Let m > be sufficiently small and assume that V fulfills 
Hvvotheses \6.1\ and \6.b\ Then lim e ^o Ey m e — E Vm . 

Proof: Since Lemma 16.91 provides upper bounds on e and we have 

lower bounds on the spectra of H Vm E which are uniform in m and e, we can 
apply Lemma 16.71 and some straightforward variational arguments to prove the 
assertion. We only point out one subtlety: In order to show that Ey.m ^ 
Ev,m,e + o(e°) we have to pick a test function in Ran(l(_ OC:A ] (Hy m £ )), for 
some A G (Ev, m ,e , e)- To ensure that this is possible and in order to have 
an e-independent bound on the numbers (Ej^ e — A) -1 entering into the constant 
in (|6.27p we need the lower bound (|6.25|) on the binding energy Y} m e — Ev,m,e c 
which does not depend on e. Without an estimate on the binding energy we 
still got Ev, m .e ^ Ev, m + °( £ °)> but we would not have a useful lower bound on 

□ 

6.6 Proofs of the existence of ground states with mass 

The next theorem asserting compactness of spectral projections of H Vm asso- 
ciated with sufficiently low energies is the final result of this section. As in [4] it 
is proved by estimating the trace of the spectral projection from above by the 
trace of some finite rank operator, namely the one in (|6.37|) . This finite rank 
operator is constructed by means of a suitable restriction of the discretized field 
energy with discrete spectrum and a harmonic oscillator potential which com- 
pactifies the electronic part of the operator. In order to sneak in the harmonic 
oscillator potential in the proofs below we exploit the exponential localization 
of low-lying spectral subspaces once more. The latter idea stems from [SHI ES] 
and replaces an argument in [4] that works only for small e and/or A. 

Theorem 6.11 Let e e R. A e (0, oo), and assume that V fulfills Hypothe- 
ses \6.1\ and 1 6. 6\ Define \ '■= l(_oo e % +m/4]^v m o) and assume that m > 
is sufficiently small. Then Tr{\} < oo. In particular, E Vm is an eigenvalue of 

bothH V,m,0 andH l, m - 
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In the proof of the preceding theorem, which is carried through separately for 
(j = PF and ft — np below, we shall employ an orthogonal splitting of J?^* into 
subspaces of discrete and fluctuating photon states, 

Here P £ is defined in (^TT0]> . The splitting jQ = J(T£ © gives rise to an 
isomorphism 

L 2 (R 3 ,C 4 )®J^ b [^] S (L 2 (R 3 ,C 4 )® J? b [jO) ® ^b[J^], (6.31) 

and we observe that the Dirac operator and the field energy decompose under 
the above isomorphism as 

D Am , s = Da*,,. ® 1 J , Et,m,e = Hf, m ,e ® ^ + I* ® H f tme . (6.32) 

Here and in the following we designate operators acting in the Fock space factors 
&h[Jfm], £ G {d,f}, by the corresponding superscript d or /. In fact, the dis- 
cretized vector potential A mj£ acts on the various n-particle sectors in &\>[JfT>\ 
by tensor-multiplying or taking scalar products with elements from (apart 
from symmetrization and a normalization constant). 

For I e {rf, /}, we denote the identity on ^[Jt^] by \ l and the projection 
onto the vacuum sector in .^[Jt^\ by P^i and write P^ e := \ l — Pqi. The 
identity on £ 2 (R 3 , C 4 ) is denoted as l el . 

The semi-relativistic Pauli-Fierz operator. 

Proof: We prove Theorem lG. Ill with jj = PF. On account of LemmaHTTJi) 

we have X H v F m ,oX > X^y,m, e X ~ o(e°)x- Using f|6.32[) and H f i m e P iV = 0, 
we then obtain 

x{Co-C- m /2}x + o(e°)x 

> X { [ \D At c \ + V + Hf^ e - E%% - m/2] ® P w } x (6.33) 
+ X { [ \D A , mE \+V + Ht„^ - E^ £ ] ® P^ } X (6.34) 
+ X { l cl ® l d ® (< m , E - ^y,™ + ~ rn/2) } X ■ (6.35) 

In view of (IB321) we have Ey F rn £ = inf a[ \D A ^ e \ + V + Hf m e ]. Therefore, 
the expression in (|6.34[) is a non-negative quadratic form. For sufficiently small 
e > 0, the expression in (|6.35[) is a non-negative quadratic form, too, because of 
Hf m P^f ^ m P^ f and Corollary 16.101 In order to treat the remaining term 
in (|6 33|) we write 

[\D Adm J+V + H?^ e ]®P nf 
= (1 ® Par) { \D Am J + V + H f , m , e } (1 ® P Q f) 

> (1 ® Pn/ ) {e |D | + e F f , m , s - const(e, m, V, e, A)} (1 ® P w ) . 
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In the second step we assumed that e > is small enough. Altogether, we arrive 
at 

X { H^ m ,o - El ¥ m - m/2 } x + o(e°) X + e X N 2 X 

> X { ( e \Do\ + £ |x| 2 + e P f d m , £ - const ) <g> P n/ } x 

> x{ [e \D \ + e\x\ 2 +eH?^ £ - const] _®P nf } x , (6.36) 

where [•••]- ^ denotes the negative part. The crucial point about the previous 
estimate is that both \D \ + |x| 2 and Hf m £ have a purely discrete spectrum as 
operators in the electron and discrete photon Hilbert spaces. Besides Pq/ has 
rank one, of course. (w e has a discrete spectrum as an operator in Jff^ = P e Jt 
because the eigenspace in J£^J corresponding to some value attained by uj e is 
finite-dimensional. Using w £ ^ m > it is then easy to see that the spectrum 
of its second quantization, Hf m £ = dT(uj e \j^d), is discrete, too.) In particular, 
we observe that 

W~. e := [e\D \+e\x\ 2 + eHl„ h£ -const] _®P nf (6.37) 

is a finite rank operator, for every sufficiently small e > 0, no matter how large 
the value of the (e, m, V, e, A)-dependent constant is. As a simple consequence 
of the exponential localization we further know that x l x | 2 X 1S bounded. Us- 
ing x { -Hy.m.o — Ey F m — m/2 } x ^ — (m/4) x we obtain the bound (o(e°) — 
m/4) Tr{x} ^ Tr{xW~ e x} > —oo from (|6.36[) . Fixing e > sufficiently small 
we conclude Tr{x} < oo. □ 
The no-pair operator. 

Proof: We prove Theorem 16 . 1 1 1 with (j = np. On account of Lemma l6.7f i) 

we again have X-ffy- P „ l0 X ^ xHy P mE x — o(e°) x- In view of (|6.32|) we have 
P^ m e = P^ d ® 1^ with P+ d := 1[ .oo)(-Da^ J an d we observe that Py P me 
decomposes under the isomorphism (|6.3ip as 

Hy V m ,e = X2®V + Pt*^®Hl mie , (6.38) 

X " == P k. c {D ^, +V + H ^e) Pl te ■ 

Writing l cl ® l d = P+ d + P Ad and P = P n , + P£ f and using ([638]) and 
H tm,e P nf = °> we obtain 

x{H v p m , -E v * m - m /2}x + o(e°)x 

> X { - (^y P m + m/2) P^J ® P n / } X (6.39) 

+ x { (A £ d - P^ p m , £ P+J ® i& } x (6.40) 

+ X { Pt ttC ® (Kfm.e " Kf™ - + //;'„, )/V } X (6.41) 

-(£^ + m/2) X {P A , ®l f } X - (6.42) 
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Next, we observe that Ey' m E — infer [X*] and proceed as in the proof of The- 
orem 16.111 with jj = PF: We omit the expression in (|6.40[) which is a non- 
negative quadratic form. For sufficiently small e > 0, the term in (|6.41[) is non- 
negative also, since H( m e P^ } > m P^ f and lim £ ^ £yf m , E = Ey P m by Corol- 
lary [00] The term in (|6.42|) . is some o(e°) on account of xl-P/Td <& V}x 

X ( P a> - P A m JX, HSIISl), and the boundedness of e F h) 12 X - (The latter fol- 
lows from the exponential localization and higher order estimates.) Putting all 
these remarks together and applying p.32[) in the last step we obtain 

x{H v p m , -E v V m -m/2} X + o(e ) X 

> X { {X d E - (E v % + m/2) P+ t ) ®P nf }x 

= X (1 ® Pat ) {f^f w , e - (Ey p m + m/2) P+ m e } (1 ® P n/ ) X 

> X Pt... { [ £ 1^° I + £ I X ! 2 + £ H f,m,e - COnst ] - ® P H/ } i^,,. X 

-£X^ m J|x| 2 ®Po/}^ m , e X, 

where the constant in the penultimate line again depends on e, m, V, e, and 
A, and the operator in the last line is bounded due to the localization of x- We 
may thus conclude as in the end of the proof of Theorem 16.111 with ft = PF. □ 

7 Infra-red bounds 

The final step in the proof of the existence of ground states is a compactness 
argument showing that a sequence of normalized ground state eigenfunctions, 
4>mj i Trij \ 0, of operators with photon masses rrij contains a strongly conver- 
gent subsequence whose limit turns out to be a ground state eigenfunction for 
the original operator. This compactness argument is explained in the subse- 
quent Section [S] As a preparation we now discuss the infra-red bounds stated 
in the following proposition. They are proved in [28 for the semi-relativistic 
Pauli-Fierz operator and in |29j for the no-pair operator starting from a suitable 
representation of a{k) cf> m . In order not to lengthen the present exposition too 
much we only outline the proof of the soft photon bound for the semi-relativistic 
Pauli-Fierz operator in Subsection 17.21 We recall the notation 

(<z(fe) </>)W (&!,...,*;„) - (n+l)^ 2 ^ n+1 Hk,k u ...,k n ), neN , 

almost everywhere, where tp = (if}( n ')%L € and a{k) Q = 0. 

Proposition 7.1 Let e € R, A 6 (0, oo), and assume that V fulfills Hypothe- 
ses \6.1\ and \6.6l . Then there is a constant, C > 0, such that, for all sufficiently 
small m > and every normalized ground state eigenvector, <j>^ m , of H Vm , we 
have the following soft photon bound, 

2 C 
\\ a { k )4L\\ < l{m<|k|<A} TTT j (7-1) 
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for almost every k — (k, A) £ R 3 x J.2, as well as the following photon derivative 
bound, 

|| a(k) 4>L - a(p) A || < C |k - p| (^L^ + , (7.2) 

for almost every k — (k, X),p — (p,n) £ R 3 x Z2 wii/i to < |k| < A and 
to < |p| < A. 

We remark that the photon derivative bound (17.21) is actually the only place in 
the whole article where the special choice of the polarization vectors (|2.4I) enters 
into the analysis. 

7.1 The gauge transformed operator 

In order to derive the infra-red bounds (|7.1[) and (|7.2[) by the method outlined 
in Subsection 1 7. 2 1 it is necessary to pass to a suitable gauge [5J[T7]. For otherwise 
we would end up with a more singular infra-red behavior of their right hand 
sides. To define an appropriate operator-valued gauge transformation [17] we 
recall that, for i,j £ {1,2,3}, the components im(x) and Am (y) of the mag- 
netic vector potential at x, y £ R 3 commute in the sense that all their spectral 
projections commute; see, e.g., Theorem X.43 of [23] • Therefore, it makes sense 
to define 

r® 3 
U := l C 4®C/ x d 3 x, U x := nV 35 ^™ <°>, x= (x 1 ,x 2 ,x 3 ) £ R 3 . 

Then the gauge transformed vector potential is given by 

A m (x) := A ro (x)-A m (0) = a f (G m ) + a(G m ) , 
where a"(G m ) = f® 3 l c * ® a'(G Xjra )d 3 x, and 

G x , m (fc) := -e 1{m<l ^- A} ( e lk ' x -l)£(fc) - (e ik ' x - l)G , m (fc) , 

27T-\/|k| 

for all x G R 3 and almost every k = (k, A) £ R 3 x Z2. Here G x m is defined in 
(|6.1I) . In fact, using [[/ , a ■ A m ] = we deduce that 

U D Am U* = D Am , U P+ m U* = P± m , U \D Am I U* = \D Am \ . 

The crucial point observed in [6] is that the transformed vector potential A m 
has a better infra-red behavior than A m in view of the estimate 

|G x , m (fc)| < |k| |x| |G , m (fc)| . (7.3) 
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In particular, infra-red divergent (for m \ 0) integrals appearing in the deriva- 
tion of the soft photon bound are avoided when we work with A m instead of 
A m . We further set 

H t := UHfU* = Hf + ^x • (o(w G„, m ) - a) \u G , m )) (7.4) 
+ 2(ux- G ,m |x- G ,m) , 
H* F m := [/ H™ m U* = \D Am \+V + H f , 
Hy P m ■= U H^ m U* = P± m (D Am + V + Hi) P± m , 

? m ~U4t, (7.5) 

so that $ m is a ground state eigenfunction of Hy m . One can easily show that, 

if the infra-red bounds (|7.1I) and (|7.2I) hold true with $ m replaced by then 
they are valid for $ m as well (with a different constant, of course). 

7.2 Soft photon bound for the semi-relativistic Pauli-Fierz 
operator 

To simplify the notation we write <f> m instead of in this subsection. A brief 
calculation yields 

< ( ( H^n ~ Ev*m ) <*(fc) Im \ a{k) 4> m ) = ( [H™ m , a(k) } 4> m \ a(k) 4> m ) 
= ( I S A m ' a ( k ) 1 D A m (t>m\a(k)(f> m ) + ( % m [ V Xm , a(k) ] (j> m | a(fc) (j> m ) 
+ {[H h a{k)]4> m \a{k)4> m ). (7.6) 

Combining the commutation relations 

[a(k),aHf)] = f(k), [o(fc),o(/)]=0, [o(fc),fl- f ] = u(A;)o(fc). 

with the formula (|7.4j) for Hi we obtain 

[iff , a(k) ] = —uj(k) a(k) + iu(k) x • Go. m (fc) . 

Furthermore, we get [Dj^ , a(k)] = —a ■ G x ,m(&) and 

[S Am ,a(k)}= f°° R A Jiy)[a(k),D x jR A Jiy)^, (7.7) 

where i?^ (iy) = (D A — iy)^ 1 . Next, we insert (|7.7|) into (|T.6[) . move the term 
containing uj(k)a(k) to the left hand side, and divide by oj(k). Furthermore, 
we introduce a weight function, -F(x) = a |x| 2 /-\/l + |x| 2 , for some small a > 0. 
Abbreviating 

:== ef D A m e ~ F = D A m + ia ■ VF > R l (*v) ■= ( D i m ~ f f 
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we obtain the following result, 

/oo 
(R A Jty){cx-G^ m e- F \k\- 1 } x 
OO 

x fl£ m (iy) D? m e F 4> m | a(k) 4> m ) ^ 

- ( S A m i a ' G x,™ e~ F Ikr 1 } e F cj> m I a(k) 4> m ) 

- iG , m (k) ■ ( (xe~ F ) e F (j> m \ a(k) cj> m ) . (7.8) 

The purpose of the exponentials e~ F introduced above is to control the factor 
|x| coming from (|7.3p . In fact, 

||a-G x , m e- F || < V2 sn V \G^ m e- F ^\ ^Oil)^' 2 t {m ^ KK} . 

X 

Using this it is easy to see that the sum of the last two expressions in (I7.8[) is 
not greater than 0(1) |k| _1 l{ m ^|k|^A} l|e F 0m|| 2 + || ct(fc) 4> rn \\ 2 /2. By virtue of 
the second resolvent identity we further get 

Since we have (iy)\\, \\R Am (iy)\\ ^ 0(1) (1 + | | ) 1 by (12T2H and 

\\RAj^y)\DAj 1/2 \\^o(m + \y\r 1/2 , 

we thus obtain 

II^L^iy^H ^ 0(l)(l + \y\)- 1 / 2 \\\D A J 1 / 2 e F j m \\. 

Therefore, the integral in (17.81) converges absolutely and we arrive at 

\\a(k) m || 2 < 0(1) Ikl^l^ikKAj || \D A J 1/2 e F m \\ 2 + 3 \\a(k) m || 2 /4. 

If a > is small enough, then || \D A \ x / 2 e F <j) m \\ is bounded uniformly in m > 0, 
due to a strengthened version of our exponential localization estimates; see 
Lemma 5.4 of [28] whose proof works for every potential V fulfilling Hypothe- 
sis 16.11 Together with the last remark in the preceding subsection this yields 
the soft photon bound. 

The above calculations illustrate the importance of the formal gauge invari- 
ance of our models. In fact, without the gauge invariance we could perform these 
calculations only for G x . m instead of G Xjm . In this case we would, however, 
end up with an upper bound on ||a(fc) (f) m \\ 2 of the form 0(1) |k|~ 3 l{ m ^|k|<gA} j 
which is not integrable near zero and, hence, is not suitable for the arguments 
presented in the following section. 
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8 Existence of ground states 



We have now collected all prerequisites to show that inf a [H v ] is an eigenvalue 
of H v . The final compactness argument is presented in the first subsection 
below. The main idea behind it is borrowed from [T7J: Namely, to show that 
{ ( l > rn j }ji 171 j \ 0: contains strongly convergent subsequences we may restrict 
our attention to finitely many Fock space sectors and to a compact subset of 
the k-space. In fact, this is a consequence of the soft photon bounds. On 
account of the exponential localization estimates it further suffices to consider 
compact subsets of the x-space. Moreover, the photon derivative bounds and 
the localization in energy lead to bounds on the (half-)derivatives of the vectors 
<j)^ mj w.r.t. x and k on compact sets and in the finitely many Fock space sectors. 
The idea proposed in [17] is to exploit such information by applying suitable 
compact embedding theorems. Essentially, we only have to replace the Rellich- 
Kondrashov theorem applied there by a suitable embedding theorem for spaces 
of functions with fractional derivatives. (In the non-relativistic case the ground 
states 4>\n possess weak derivatives with respect to the electron coordinates, 
whereas in our case we only have Inequality (|8.5I) below. For a variant of the 
argument that avoids the Nikol'skh spaces introduced below by switching the 
roles of the electronic position and momentum spaces see [25] .) 

In the second subsection we discuss the degeneracy of the ground state en- 
ergies by applying Kramers' degeneracy theorem similarly as in |40) . 

In the whole section the coupling function G x is the physical one defined in 
(& 

8.1 Ground states without photon mass 

We shall apply the following elementary fact: 

Let S be some non- negative operator acting in some Hilbert space, X . Let 
{r]j}j£M be some sequence in , converging weakly to some n £ such that 
Vj € Q(S) and (rjj \Srjj) = US' 1 / 2 rjj\\ -> 0. Then r\ belongs to the domain of S 
and S i] = 0. 

In fact, the linear functional /(</>) = (t) 1 S 1 / 2 4>), 4> G V^S 1 / 2 ), satisfies 
/(<£)= lim(7 7 ,|S' 1 / 2 0) = lim(S 1 / 2 % |0)=O, 

j— too j— too 

and the self-adjointness of S 1 / 2 implies n £ V(S 1 ^ 2 ) and S 1 / 2 r) = 0. 
We remind the reader of the notation (j £ {np, PF}. 

Theorem 8.1 Let e £ R, A > 0, and assume that V fulfills Hupotheses \ 6.1\ 
and \6.(h Then E v is an eigenvalue of H v . 

Proof: For every sufficiently small m > 0, there is some normalized ground 
state eigenfunction, 4>lm 01 H Vm . We thus find some sequence nij \ such that 
WmjIjetN converges weakly to some cjfi £ . According to Theorem 13.41 the 
form domains of Hy F and Hy m , m > 0, coincide whence £ Q(Hy F ). By 
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the characterization of the form domain of the no-pair operator in Theorem l3.6l 
we further know that P£ (f)^. G Q(Hy P ). It is trivial that Pj^^m- converges 

weakly to Pt (f> np in F+ Jt. Since Nj := (Pt - P A ) H~ 1/2 converges to zero 
in norm, where Hf = Hi + E, for some sufficiently large E 1, we further 
have ( ip | P A np > = linij^co( Nj h\ 12 | ) = 0, for every i> G 0, whence 
PT </> np = 0, or, PT np = np . Hence, by the remark preceding the statement 
it suffices to show that $ ^ and ( ^ | (if*. - J5*.) ^ ) -> 0. 

The latter condition immediately follows from lim™^ Ev,m = Ev and 
Lemma l6.5f ii) which together imply 

( fa. | (Hi - E v ) ^ ) = ( ^ | (< mj - < mj ) ^ ) + o(m°) = o(m°) . 

In order to verify that f / we adapt the ideas of [T7] sketched in the first 
paragraph of this section. 

We write <p m = (<t> { 2\)n=o G 0^=0^ P*1 in what follows - Let £ > 0. By 
virtue of the soft photon bound we find uq G N and C G (0,oo) such that, for 
all sufficiently small m > 0, 

E ii^Sn 2 < ^E-ii^an 2 = -f- / iK*)dLii a « < f < f ■ (8.i) 

* — ' Tin — Tin / Tin ^ 

n=n U n=0 U J U 

By the exponential localization estimates of Theorem 15. 1[ which hold uniformly 
for small m > 0, we further find some R > such that 

/ ll^mllcW*)* 3 * < I' ( 8 - 2 ) 

In addition, the soft photon bound ensures that <^l(x, fci, . . . , k n ) = 0, for 
almost every (x, q, k u . . . , k n ) G R 3 x {1,2,3,4} x (R 3 x Z 2 ) n , k 3 = (k,-,Aj), 
such that |kj| > A, for some j G {1, . . . ,n}. (Here and henceforth <; labels the 
four spinor components.) For < n < uq and some fixed 9 — (?, Ai, . . . , A„) G 
{1,2,3,4} x Z£ we set 

^m,?,tt ( ' X ' kl ''"' kn ' ) " ^m,tt( X ' ? ' k l' A l'---' k n' A ™) 

and similarly for 0' = (<M )5£Ln- Moreover, we set, for every (5^0, 

Q n< s := {(x,ki,...,k n ) : |x| < i? - 5 , 6 < \kj\ < A - 6 , j = 1, . . . , n } . 

Fixing some small 6 > we pick some cut-off function X € Cg°(IR 3 ( n+1 ), [0, 1]) 
such that x = 1 on Q n ,2S and supp(x) C Q n .8 and define ^ := X «• ^ s 
a next step the photon derivative bound is used to show that {V* 37 .}me(oi] 
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is a bounded family in the anisotropic Nikol'skh space@ i^(R 3(n+1) ), where 
s = (1/2,1/2, 1/2,1,..., 1) and q = (2, 2, 2,p, . . . ,p) with p € [1,2). In fact, 
employing the Holder inequality (w.r.t. d 3 xd 3 ( n_:L )K) and the photon derivative 
bound (|7T2"1) . we obtain as in [T7], for p G [1, 2), m € (0, 8], and h € R 3 , 

0j» if (x, k + h, K) - ^(x, k, K)|" d^kd^K 

{5<|k+h|<A} 

< C X] / || «(k + h, A) 0» n - a(k, A) 0»„ fd 3 k < C |h|", 

-^6 Z 2 m <|k|<A, 

771<|k+h|<A 

where the constants C,C G (0,oo) do not depend on m G (0,(5]. Since </>^^ 
is symmetric in the photon variables the previous estimate implies |42[ §4.8] 

(n) 

that the weak first order partial derivatives of <b - with respect to its last 3n 
variables exist on Q n< s and that 

n 3 

i=i »=1 

for m € (0, 5] and some m-independent C" € (0, oo), with r := (0, 0, 0, 1, ... , 1). 
The previous estimate implies HV^^ J W£(R 3 <™+ 1 >) ^ C"> f° r some C" € (0,oo) 
which does not depend on m G (0,5]. Moreover, the anisotropic Sobolev space 
W p r (R 3 (™ +1 )) is continuously embedded into fl^R 3 ^ 1 )); see, e.g., g2J §6.2]. 
By Theorems 13.41 and 13.61 we get, for rt G N, 

c _1 < <fi I I A> | <fi ) < ( 4> m ,i I < m ^,8 > + c = < m + c < 4 + 2c , (8.5) 

for some m-independent c G (0, oo). Therefore, {vj^jj Am&(p,S\ anci : hence, 
{V'^g j}me(o ,5] are bounded families in the Bessel potential, or, Liouville space 
L|'(R 3 ( n+1 )), r / := (1/2,1/2,1/2, 0,...,0), where the fractional derivatives are 

2 For n,...,r d 6 [0,1], qi,...,q d > 1, wc have tf&.7.'/<£ <,) (R d ) := nti^,'x,(R' ! ). For 
r,: 6 [0, 1), a measurable function / : R d — > C belongs to the class . (R d ), if / e L^R" 1 ) 
and there is some M" g (0, oo) such that 

\\f(-+hei)-f\\ L « i(Kd) H M\hp, hen, (8.3) 

where is the i-th canonical unit vector in R d . If rt = 1, then 118.31 1 is replaced by 

\\f(- + hei)-2f + f(--hei)\\ Ln(Rd) ^M\h\, heR. (8.4) 

Hq^'l^q^ (R d ) is a Banach space with norm 

ll/lfcfi'* : = ^ll/llrwCR*) + ^<!? d M *' 

where Mi is the infimum of all constants M > satisfying 118.31 1 or 1 18. 4> . respectively. Finally, 
we abbreviate H q T1 '-' Td) (R d ) := H^. 1 .:;^ d) (R d ). 
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defined by means of the Fourier transform. The embedding L\ (R 3 (" +1 )) — > 
#r'( R 3(rs+i)) is continuous, too; see §9.3 in [42]. Altogether it follows that 
j}me(o,5] is a bounded family in i/q(IR 3 (™ +1 )). Now, we may apply The- 
orem 3.2 in [IT]. The latter ensures that Ame(0,5] contains a sequence 
which is strongly convergent in L 2 (Q n , 2 s) provided 1 — 3n(p^ 1 — 2 -1 ) > 0. 
Of course, we can choose p < 2 large enough such that the latter condition 
is fulfilled, for all n = 1, ...,no — 1. By finitely many repeated selections of 
subsequences we may hence assume without loss of generality that {</>|™^ -g ^Ijefti 

converges strongly in L 2 (Q nt2 g) to its weak limit <Pgli f° r ^ n < no. In 
particular, by the choice of n and R in (|8.1I) and (|8.2I) . 

ii^ii 2 > A- EE > A- IK ii 2 - - °^°) . 

n=0 e 

where we use the soft photon bound to estimate 

n -l 

E ElkLt^ 1 ^ n<2*v a- 25} || 2 

n=l g 

< £ | ||a(k,A)^.H 2 rf 3 k = o(5°), 6\0. 

■^£ Z 2 {|k|sJ25}U 
{|k|^A-23} 

Hence, ||0"|| 2 ^ 1 — e — o(6°), where 5 > and e > are arbitrary, that is, 
1 10" 1 1 = 1. (In particular, $ m . $ strongly in □ 

8.2 Ground state degeneracy 

Suppose that V(x) = V(— x). As already mentioned in the introduction it is 
remarked in [40] that every (speculative) eigenvalue of Hy F and, in particular, 
its ground state energy is evenly degenerate in this case. The authors prove this 
statement by constructing some anti-linear involution commuting with Hy F and 
applying Kramers' degeneracy theorem. We shall do the same for the no-pair 
operator in the next theorem which originates from [29) . 

Theorem 8.2 Let e € R, A > 0, assume that V fulfills Huvotheses lfi.il and 
assume in addition that V(x) = V(— x), for almost every x. If the ground state 
energy E v p is an eigenvalue o/iJy P , then it is evenly degenerate. 

Proof: Similarly as in [40] we introduce the anti-linear operator 

d := Ja 2 CR = —a 2 JCR, J := 

where C : Jf — > denotes complex conjugation, C -0 := i> , i> & and 
R : Jff — > Jff is the parity transformation (Rip) (it) := tp(— x), for almost every 
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x G R 3 and every V € JT = £ 2 (R|,C 4 ® Obviously, [<!?, -iS^] = 

[$ ,V] = [i? , iJf] = 0, on P(|Do|) n V(Hf). Since «2 squares to one and C«2 = 
— «2 C, for all entries of a 2 are purely imaginary, we further get {f 2 = — 1 and 
02] = 0. Moreover, the Dirac matrices ao, ot\, and 0:3 have real entries 
and [JQ2 , ctj] — J {a 2 , ctj} — by ()2.1[) . whence [# , ccj] = 0, for j e {0, 1, 3}. 
Finally [1?, e ±ikx ] = implies [1?, A«] = on V{H] /2 ), for j G {1,2,3}. 
It follows that [d . Da] = on f = i)f and, since Da is essentially self- 
adjoint on 9, we obtain tfV(D A ) = D(Ka) and [1?, D A ] = on V(D A ), 
which implies $ RA.(iy) — Ra{—w)^ = on Jif, for every y £ R. Using the 
representation (|2. 18[) we conclude that [d , P£] = on JO?. In particular, $ can 
be considered as an operator acting on J^a ■ Furthermore, we obtain H v p 'dip — 
i9 Hy P tp = 0, for every jj£f. Since ^ is a form core for ii^ p we can easily 
extend this commutation relation and show that 1? maps T>(H V P ) into itself and 
Hy P, dip = ■dH v p tfj, for every ip G V(H V P ). Hence, by Kramers' degeneracy 
theorem every eigenvalue of H v p is evenly degenerated. (In fact, H v p <f> — Ey <f> 
implies H v p tf(f> = E v •&</>, and 0_Ltf0, since (tf(f>\(f>) = -{ 1? 1 1? (1? 0) ) = 

-<c>|e^) = -<<?</> I </.).) □ 

Using similar arguments we derive the next lemma which has been referred to 
in Sections [5] and [6l We use the notation introduced in Subsection 16 . 1 1 and write 

for the analogs of the no-pair operators on the negative spectral subspace which 
appeared in (|5.ip and (|6.28|) . Then H' ip me = HQ P me + i?Q^~; see (|5.ip and 
(|6.28p . We unify the notation by setting Hq q := Hq P , etc. 

Lemma 8.3 Let e G IR, A G (0, 00), and m > e 0. Then 

SSf„ d = f - infal^J = infa^-] = infafej. 

Proof: Let C and i? be defined as in the preceding proof. We introduce the 
anti-linear operator r : ~~ ^ ^m"; r := ot^C R. Similar as in the previous 
proof we verify that tDa„ 1e = --Da„ J on T>(Dp^ m E ), thus RA me {iy)T — 
— t RA m , s (iy), thus P Am £ T = t Pjf again by (|2.18|) . Notice that we use the 
property (|6.9p of the discretized phase in A mj£ to obtain these relations in the 
case e > 0. Consequently, Hq p m e r = r ffj^~ on a natural dense domain (again 
called in with r£F = By Theorem 13.71 we may assume that the no- 
pair operators in the plus or minus spaces are essentially self-adjoint on Pj^ m @, 
respectively, and we readily conclude. □ 

9 Commutator estimates 

In this final section we collect some bounds proved in |37j on the operator 
norm of various commutators which have been used repeatedly in the preceding 
sections. In the whole section we assume that G x fulfills Hypothesis 12. II 
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9.1 Basic estimates 

Our estimates on commutators involving the field energy Hf are based on the 
next lemma. The following quantity appears in its statement and in various 
estimates below, 

Si = 6 U (E) 2 := 8 / Wv{k ;^ ? ||G(fc)||L dk, E,v > , (9.1) 

where w v (k, E) := E 1 ' 2 ^ {{E+us{k)) v+1 / 2 -E v (£+w(fc)) 1/2 )- We observe that 
w 1 / 2 (k,E) < uj(k) and, hence, Si/ 2 (E) < 2di, E > 0. Moreover, < &/(l), 

for _E 1. We recall the identity 

< Hl /2 cj) | i? f 1/2 V ) = / w(k) ( o(A)^ | a(k)tp )dk , 0, V G V(H\ j2 ) , (9.2) 



which is a consequence of the permutation symmetry and Fubini's theorem. 

Lemma 9.1 Let v, E > 0, and set Hi := Hi + E. Then 

|| [a ■ A , Hf v ] H'( || < S V (E)IE X I 2 . (9.3) 

PROOF: ([37]) We pick <f>, tp e and write 

(^[a-A.tff]^) 

= (^|[a-a(G),^]^V> - ([a-aCG),^]^!^^)- (9-4) 

By definition of a(fc) and iJf we have the pull-through formula a(k) 6(Hf) ip = 
9(H{ + ui(k)) a{k) ip, for almost every k and every Borel function 9 on R, which 
leads to 

[a{k),HY v ]H^ 

= { ((H { + u(k))- u - Hf") (Hf + u{k)Y +1 ' 2 } a{k) H~ 1/2 V . 

We denote the operator {• • • } by F(k). Then F(k) is bounded and 

d (t + E + bj{k)) u+1 / 2 , 
as 



\\F(k)\\ < / sup 
Jo t^o 



ds (t + E + su(k)) u 



l 1 ^^ / \ s = Mk,E) / E^. 
J ds (E + su(k)Y K " 

Using these remarks, the Cauchy-Schwarz inequality, and (|9.2I) . we obtain 
\(4>\{a-a(G), Hf v ]H^)\ 

< J \\<j>\\ \\a ■ G(k)\\ \\F(k)\\ \\ a(k) Hf 1 ' 2 \\ dk 

< 11011 (2 |^^||G(fc)||Ldfc) 1/2 ( J u(k)\\a{k)H; 1/2 ^\\ 2 dk) 1/2 

$v(E) .. ., I, 1/2 £—1/2 , II 
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A similar argument applied to the second term in (|9.4j) yields 

|( [a • a(G) , Hi") <j> \ H? i> > | < ^§ || ff f 1/2 # f " 1/2 || ||V>|| , 

where 8„(E) is defined by (|9 . 1 j) with w u (k,E) replaced by 

w u (k, E) := E 1/2 - v {E" (E + tu(k)) 1/2 - E 2v (E + lu) 1/2 -") . 

Evidently, w v ^ w v , thus 5 V ^ S u , which concludes the proof. □ 
Choosing E large enough, we can certainly make to norm in (I9.3P as small as 
we please. This observation can be exploited to ensure that certain Neumann 
series converge in the next proof which yields convenient formulas allowing to 
interchange the field energy with resolvents of the Dirac operator. 

We define J : [0,1) -> R by J(0) := 1 and J(a) := V6/(l-a 2 ), for a e (0,1), 
so that \\R F (iy)\\ < J(a) (1 + y 2 )~ 1/2 , where R° A (iy) := R A (iy); recall (j^20|) 
and (p^Tj) . 

Corollary 9.2 Let y e R and F G C°°(R^, R) sucft tfiai |VF| < a < 1. Assume 
i/iai z^, E > satisfy 5 V J{a) / E 1 / 2 < 1, and introduce the following operators, 

T v := [fff , a ■ A] fl?, 

oo oo 

(%) := T,F, T? (iy) := £{-2? flI(tj/)F". (9.5) 

J=0 ]=0 

Then \\T U \\ ^ 8,/E 1 ' 2 , \\Q F {iy)\\ s? (1 - 6 V J{a) / E 1 ' 2 )' 1 , 6 G {S,T} ; and 

i?r ^a(^) = (iw) Rl(iy) 8f" , (9-6) 

«A(*y) Hr = Hr Rliiv) Tf (itf) , (9.7) 

In particular, R A (iy) maps T>(Hf) into itself. 

Proof: The following somewhat formal computations show the simple 
idea behind the proof (see Section 3 of [37] for more details), 

Hf v R F A (iy) = R F A iiy) Sf + R F (iy) [D A + ia-VF-iy, Hf] R F (iy) 
= R F A (iy) Hf" + Rliiy) [a ■ A, Hf] Hf(Hf u flftiy)), 

which implies that R^(iy) Hf" = (l + R A (iy) T„) _1 Hf u R^(iy). The oper- 
ator inverse appearing here is given by the Neumann series "E. F {iy) when we 
choose E in Hi = H s + E so large that 5 v J(a)/E 1 / 2 < 1. □ 

9.2 Commuting projections with the field energy 

Lemma 9.3 Let a, k G [0,1), let v, E, and F be as in Corollary \9. 6 A and 
assume that F is bounded. Define C F := e F [S Al Hf"]Hfe~ F on T)(Hf). 
Then 



58 



In particular, Sa maps 'D(Hf) into itself and, if E x l 2 > 5 U) then the following 
identities hold true on T>(H\), where C v := (C°)* £ JtC(Ji?), 

Hf Sa — Sa H{ + C u H\ , 5a Hf = H\ Sa + H\ C* . (9.9) 

Proof: ([37]) Combining flM]) with (|2~Tgj) we obtain, for all <j) £ V{\D A \ K ) 
and ijj £ V{H%), 

\(\D A \ K cb\C F ^)\ < / | < |i> A r ^ I r„ iZi(<v) ^ > | — 

< ||T„|| / ||pAri2i(i!/)0||||Sf(<i/)|||K(iy)||^|^||W. 

JtR K 

Here we estimate ||T„|j by means of (|9.3I) and we write 

\D A \ K Rl(iy) = \D A \ K R A (iy)(l-ia-VFRi(iy)), 

where || |-D A | re E a (m/)|| «S const(«)(l + y 2 )- 1 /^/ 2 . Moreover, ||Sf(iy)|| sC 
(1 — 5 v J{a)/E 1 / 2 )^ 1 , y £ R, by Corollarv l9.2l Altogether these remarks yield the 
asserted estimate. Now, the identity Sa i?f _I/ = Sa — Hf v (C°)* in Jzf(^) 
shows that Sa maps the domain of H" into itself and that the first identity in 
dHJ) is valid. Taking the adjoint of (JOJ) and using [iT^ S A }* = S A Hf (which 
is true since Sa is densely defined and Sa = S^ 1 £ Jz? (Jf )) we also obtain 
the second identity in (|9.9p . □ 

We can now prove an estimate asserted in the proof of Theorem 13.41 namely 
T = Re6 < e|-D A | + const d 1 /E 1 / 2 , for e e (0,1] and E ^ (4di) 2 , where 
: = [|-D A |,.ff f 1 ^ 2 ]i? f 1 ^ 2 . In fact, this follows easily from Corollary 19.21 and 
Lemma E3] since 9 = |£> A | 1/a SA{|£A| 1/a C? /2 } + T 1/2 C° 1/2 - T 1/2 S A . 

9.3 Double commutators 

Lemma 9.4 Let a £ [0, 1) and let F satisfy (|5.8p . Moreover, let v, E > such 
that 8 V J(a)/£ 1/2 1/2. T/ien 

|| \D A \ [xi e F , [P£ , X2 e- F ]] || < J(o) J] ( a + Halloo) , (9.10) 

i=l,2 

||^[xie F , [P+, X2 e- F ]]i? f -'1 <8J(a) n ( a + H V ^lloo) - (9.11) 

i=l,2 

Moreover, for every self-adjoint multiplication operator V in L 2 (R 3 ) satisfying 
H 1 (R :i ) C T>(V), we find some constant CV £ (0, oo), depending only on V, 
such that 

\\\V\[xie F ,[P+,X2e- F }]Hi 1/2 \\^C v J(a) J[ (a + || V Xl ||oo) . (9.12) 

z=l,2 

In (|9.12l) we also assume that E ^ (4d\ J{a)) 2 and E ^ 1. 
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Proof: 



7\) Let 4>, ip € 3>, \\4>\\ = H^ll = 1- First, we derive a bound on 



>-<jl,Tp 



\D f 



H? [ X ie F , [R A (iy) , X2 e~ F ] } &f </> 



2tt 



Expanding the double commutator we get 

[xie F , [iJA(iy)> X2e~ F ]] = vixi, X2, F ;y) + v(X2,Xi,~F ;y) , 

where 
»7(Xi,X2,F;?/) 

:= Aa(%) a • (Vxi + Xi VF) e F i? A (*2/) a • (V X 2 - X2 VF) i? A (*2/) • 



We obtain 



D A \ R A (iy) r° v (iy) a ■ (V X i + xi VF) x 
x RF(iy) T F (iy) a ■ (V X2 - X2 VF) R A (iy) T° v (iy) j, 



2tt 



( fl + ||Vxi||)(a+||Vx2| 

(l-^/FV2)2 



J(a) 



l-5 v J{a)/E^ J v 27r(l + 2/ 2 ) 



dy 



(9.13) 



A bound analogous to (|9 . 1 3[) holds true when the roles of xi and %2 are inter- 
changed and F is replaced by —F. Consequently, 1^ is bounded by two times 
the right hand side of ([9T3)) . Altogether this shows that (j9~10| and (l9TTj) hold 
true. (Just ignore \D\\ or H{, respectively, in the above argument.) By the 
closed graph theorem we further have || |F|V>|| 2 < const(V) (||V-0|| 2 + H^ll 2 ); 
ip E iF^R 3 ). Hence, (|9~T2]l follows from (|QU]) . (|9TT|) . and the inequality 

|| \V\ipf < const(V)(|| |£> A |^|| 2 + ll^f^^ll 2 ) , <p e 9, 

which holds true, for E ^ d\, by ([321), (|3~23|) . and 

|||^M| 2 ^const(U)(||Vb]|| 2 + |M| 2 ), 



□ 
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